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Part  I:  Tectonics  and  Wave  Propagation  jn  Central  Asia 


TECHNICAL  SUMMARY 


1.  Surface-wave  analysis  proves  to  be  an  effective  way  of  determining 
focal  mechanisms  for  earthquakes  with  relatively  small  magnitudes.  By 
comparing  the  surface  waves  of  one  event  to  another  nearby  event,  the 
effects  of  the  epicenter-to-station  path  can  be  largely  reduced.  This  is 
true  even  when  most  paths  cross  or  lie  along  zones  of  complicated  defor¬ 
mation.  Two  events  in  the  northern  Arabian  Sea  have  similar  mechanisms  and 
can  be  interpreted  as  earthquakes  along  transform  features  between  the 
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Indian  and  Arabian  plates. 


A.  Focal  Mechanisms  of  Small-Magnitude  Earthquakes  in  Central  Asia 
Determined  from  Surface-Wave  Analysis 


Determination  of  focal  mechanisms  for  small-magnitude  earthquakes  in 
Central  Asia  aids  in  the  understanding  of  the  tectonic  framework  that 
characterizes  this  area.  For  earthquakes  with  magnitudes  greater  than 
about  5.5,  a  plot  of  first-motion  data  on  a  focal  sphere  usually  provides 
adequate  constraint  on  the  nodal  planes.  For  smaller  magnitude  events, 
however,  first-motions  are  not  well-recorded  at  a  sufficient  number  of 
stations  to  gvve  a  well-determined  solution.  Surface  waves,  on  the  other 
hand,  are  usually  well  recorded  to  distances  of  approximately  70s  for 
events  with  magnitudes  as  low  as  A. 5.  As  the  surface  waves  contain  infor¬ 
mation  on  the  depth  and  source  parameters,  it  is  thus  possible  by  proper 
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analysis  to  determine  the  focal  mechanisms  and  depths  of  earthquakes  with 
relatively  small  magnitudes. 

Under  this  contract  an  analysis  was  made  of  the  focal  mechanisms  of 
two  earthquakes  located  in  the  northern  Arabian  Sea.  One  event  occurred  at 
the  northern  end  of  the  Owen  fracture  zone  (30  March  1966,  *  5.3);  the 
other  occurred  along  the  Murray  ridge  about  300  km  to  the  east-northeast  of 
the  first  event  (9  November  1968,  Mg  =  5.1).  These  shocks  are  located  in  a 
region  of  complicated  tectonics  near  the  triple  junction  for  the  Indian, 
Arabian,  and  Eurasian  plates.  The  results  of  the  surface-wave  analysis 
show  that  the  nodal  planes  for  the  two  events  are  similar. 

Method  of  Analysis.  The  method  of  analysis  used  in  this  study  is  a 
comparative  event  technique.  The  observed  difference  in  phase  and  the 
ratio  of  amplitude  spectra  between  two  nearby  events  are  compared  to 
theoretical  values  to  determine  simultaneously  the  strike,  slip,  dip, 
depth  and  moment  of  both  events.  In  addition,  the  amplitude  spectrum  of 
each  event  should  agree  with  the  respective  theoretical  spectrum,  and  the 
final  solution  should  not  grossly  violate  whatever  limited  first-motion 
data  are  available.  In  this  method,  which  has  been  most  extensively 
developed  by  D.  Weidner  and  A.  Kafka  (Weidner,  1972;  Kafka  and  Weidner, 
1979;  Kafka,  1980),  the  effects  of  epicenter-to-station  paths  on  the 
observed  spectra  are  largely  eliminated  since  the  paths  are  nearly 
identical.  If  the  focal  mechanism  of  one  of  the  two  events  being  compared 
is  independently  known,  this  event  in  essence  provides  path  corrections 
for  the  second  event,  and  the  method  can  be  termed  a  master-event  tech¬ 
nique. 

For  each  earthquake  copies  of  seismograms  were  made  from  WWSSN  film 
chips  for  stations  with  recorded  surface  waves.  An  attempt  was  made  to 
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obtain  records  from  the  same  station  for  both  shocks  in  as  many  cases  as 
possible.  The  seismograms  were  then  digitized  and  interpolated  at  equal 
intervals  of  slightly  less  than  1  sec.  Next  the  group  velocity  dispersion 
was  determined  using  the  technique  of  Dziewonski  e£  al_.  (1969)  and  this 
information  was  used  to  filter  the  digitized  record  with  a  filter  designed 
to  pass  those  frequencies  arriving  at  a  particular  time  (i.e.,  time- 
variable  filtering,  Landisman  et  al.,  1969).  The  values  of  amplitude  and 
phase  at  5  frequencies  were  obtained  from  the  fourier  transform  of  the 
filtered  record.  The  values  of  amplitude  are  corrected  for  instrument 
response  and  attenuation  (V  =  1.5  x  10**);  while  values  of  phase  are 
corrected  for  instrument  response,  propagation,  and  the  origin  time  - 

digitization  time  delay. 

To  compare  the  observed  values  to  theoretical  values  three  parameters 
are  defined: 

rw  (A*  obs  -  A*  V 

o  2  =  — - £ - £ - 

<J>  EW± 

EW.(log  A°bS  -  log  AC)2 

o  2  *  - 

A  EWt 

2  EW  (log  (A  °bS/A  °bS)  -  log  (A  C/A  C))2 

n  —  -  ■  - —  ^  - 
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Superscripts  obs  and  t  refer  to  observed  and  theoretical  values, 
respectively;  Aj  refers  to  the  amplitude  of  the  jth  event;  Ao  is  the 
differential  focal  phase;  and  the  W's  are  weights.  The  summations  are 
taken  over  the  appropriate  stations  and  periods.  The  weights  are  taken  as 
proportional  to  the  theoretical  values  of  amplitude. 

When  the  theoretical  and  observed  values  are  close,  the  residual 
parameters  defined  above  are  small.  The  best  solution  is  found  by 
searching  through  the  8-dimensional  space  defined  by  the  strike,  slip, 
dip,  and  depth  of  both  events  to  find  that  set  of  source  parameters  that 
minimizes  the  residual  parameters. 

Results  for  the  Northern  Arabian  Sea.  The  two  earthquakes  that  form 
the  event-pair  for  the  northern  Arabian  Sea  are  associated  with  the  Owen 
fracture  zone  and  the  Murray  ridge,  respectively  (Figure  1).  Pertinent 
information  on  these  events  is  listed  in  Table  1.  The  azimuthal  distri¬ 
bution  of  stations  at  which  records  were  obtained  for  these  events  is  shown 
in  Figure  2.  The  stations  are  well  distributed  except  to  the  south  in  the 
Indian  Ocean.  Our  theoretical  calculations  of  Rayleigh  wave  radiation 
patterns  are  based  on  the  work  of  Saito  (  1967),  and  we  have  used  the 
oceanic  earth  model  of  Harkrider  and  Anderson  (1966)  for  the  source  region 
in  all  cases.  Several  authors  (e.g.  Weidner,  1972;  Patton,  1978)  have 
found  that  variations  in  the  assumed  earth  model  do  not  seriously  affect 
the  Rayleigh  wave  radiation. 

Event  1  has  been  previously  studied  by  Sykes  (1970).  He  determined  a 
fault-plane  solution  based  on  first-motion  readings  and  the  polarization 
of  S-waves.  His  solution,  although  classified  as  of  poor  quality,  is 
consistent  with  predominantly  strike-slip  motion  in  a  right-lateral  sense 
along  the  Owen  fracture  zone.  Figure  3a  presents  the  limited  first-motion 


data  for  this  event. 


Figure  1.  Tectonic  setting  for  the  two  earthquakes  studied.  The  locations  of  the 
events  are  indicated  by  the  triangular  symbols:  1-  30  March  1966,  2-  9  November 
1968. 


Q-30  MAR  66  X-9N0V68 


Figure  2.  Locations  of  the  stations  for  which  seismograms  were  studied  The 
paths  aiong  which  surface  waves  travelled  to  each  station  are  shown  by  the 

c“ed  ^  that  m“y  PathS  Cr°SS  °P  4l°ng  ZOnes  °f  compli* 


Figure  3.  (a.)  and  (b.)  Observed  first-motions  for  the  events  studied.  Solid  sym¬ 
bols  represent  compressional  first-motions  while  open  symbols  represent  those 
with  a  dilatational  character;  crosses  indicate  nodal  arrivals.  Smaller  symbols 
denote  data  of  poorer  quality.  A  lower-hemisphere  projection  is  used.  • 
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The  results  fcr  event  1  based  on  the  surface  wave  analysis  agree 
fairly  well  with  Syk.es'  solution.  Figure  da  shows  the  values  of  depth, 
strike,  slip,  and  dip  that  m-inimize  the  three  residual  parameters.  Most  of 
the  minima  are  well  defined  and  occur  at  similar  values  for  all  three 
residual  parameters.  The  phase  information  allows  one  to  resolve  the 
ambiguities  in  slip  and  dip  that  occur  when  only  amplitude  information  is 
considered.  In  Figure  5  the  observed  amplitude  as  a  function  of  azimuth  at 
a  period  of  40  sec  is  compared  to  theoretical  values.  The  agreement  is 
good.  Results  for  event  1  are  summarized  in  Table  2. 

The  second  event  occurred  along  the  Murray  ridge.  The  Murray  ridge  is 
a  bathymetric  feature  extending  in  a  northeast  direction  from  the  northern 
end  of  the  Owen  fracture  zone  to  the  continental  shelf  of  the  Indian 
subcontinent.  It  is  located  in  a  region  of  complicated  tectonics  near  the 
triple  junction  for  the  Arabian,  Indian,  and  the  Iranian  portion  of  the 
Eurasian  plates.  Although  the  tectonic  significance  of  the  Murray  ridge  is 
not  known,  McKenzie  and  Sclater  (1971)  have  speculated  that  there  may  be  a 
component  of  spreading  associated  with  the  feature. 

The  results  of  the  surface  wave  analysis  for  event  2  are  shown  in 
Figure  4b.  The  orientations  of  the  nodal  planes  are  similar  to  those  for 
event  1.  In  Figure  6  the  observed  and  theoretical  values  of  amplitude  are 
shown  as  a  function  of  azimuth.  As  for  the  first  event,  the  agreement  is 
quite  satisfactory. 

In  Figure  7  the  observed  differential  phase  between  tae  two  evt  :  ■=  is 
compared  to  theoretical  values  for  the  solutions  derived  from  the  surface 
wave  data.  The  observed  data  matches  well  the  theoretical  values, 
especially  considering  that  these  events  are  separated  by  a  distance 
v  "v  350  km)  that  is  large  for  the  comparative-event  style  of  analysis.  A 


(A,/A2 


oincide  quite  as  well  for  this  shock  indicating  a  solution  of  some- 
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Figure  5.  Comparison  of  the  observed  amplitude  data  with  theoretical  values  for 
the  final  solution:  30  March  1968,  4-0  sec.  The  fit  between  observed  and  theoreti¬ 
cal  values  is  quite  good.  A  better  fit  might  be  possible  if  more  weight  were  given 
to  the  longer  periods  (40,  50  sec)  at  the  expense  of  the  shorter  periods  (20.  25, 
30  sec)  in  the  computation  or  residual  parameters.  The  longer  periods  appear 
to  be  more  stable  and  contain  less  noise. 
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Figure  6..Comparison  of  the  observed  amplitude  data  with  theoretical  values  for 
the  final  solution:  9  November  1988,  40  sec.  The  fit  is  satisfactory  although  not 
as  good  as  for  the  earthquake  of  30  March  1966. 


FREQ  (SEC'1) 


Figure  7.  Observed  and  theoretical  values  of  phase  difference  for  the  event-pair 
studied.  At  low  frequencies  (<0.3  sec'1)  the  agreement  between  the  two  values 
is  good.  Phase  difference  values  for  five  other  stations  are  not  stable  and  are 
not  shown. 
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number  of  stations  are  not  shown,  however,  as  the  scatter  in  the  observed 
dif lerential  phase  inJicated  that  this  parameter  was  not  well  resolved  at 
those  stations. 

The  observed  values  of  amplitude  ratio  are  compared  in  Figure  8  to 
those  determined  theoretically.  The  theoretical  ratio  is  not  widely 
varying,  except  at  nodes,  since  the  mechanisms  are  similar.  It  is 
difficult  to  determine  visually  the  quality  of  the  fit  to  the  observed 
values,  but  this  theoretical  ratio  does  minimize  the  residual  parameter 
for  amplitude  ratio  (Figures  5  and  6). 

Conclus ions .  The  comparative-event  tecnnique  for  the  use  of  surface 
waves  to  determine  source  parameters  works  adequately  even  when  many  of  the 
epicenter-to-station  paths  cross  regions  of  complicated  tectonics.  Almost 
all  of  the  paths  in  this  study  either  cross  or  fall  along  the  Alpide- 
Himalayan  belt  of  deformation.  This  is  in  contrast  to  the  predominantly 
oceanic  character  of  most  paths  used  in  previous  studies. 

The  nature  of  the  Murray  ridge  is  not  revealed  on  the  basis  of  the  one 
focal  mechanism  for  this  region.  Although  a  mechanism  indicating  active 
spreading  was  not  found,  this  possibility  cannot  be  ruled  out.  The  earth¬ 
quake  studied  may  represent  movement  on  a  transform  fault  connecting  small 
segments  characterized  by  spreading.  In  this  interpretation  the  north¬ 
easterly  trending  nodal  plane  is  considered  the  fault  plane. 

An  alternative  explanation  involving  left-lateral  movement  on  the 
northwesterly  trending  nodal  plane  may  also  be  possible.  There  are  several 
short  scarps  (Figure  1)  that  have  a  strike  similar  to  that  of  this  nodal 
plane.  The  overall  tectonics  of  the  region,  however,  seems  to  favor  the 
first  interpretation. 


Figure  8.  Comparison  of  observed  and  theoretical  values  of  ratio  amplitude  for 
the  event-pair  studied.  Note  that  the  observed  values  scatter  within  a  small 
range  about  a  value  of  1.0.  This  signifies  that  the  two  events  have  similar 
mechanisms.  Any  set  of  similar  mechanisms  will  fit  the  data  fairly  well;  thus, 
when  the  mechanisms  are  similar  the  ratio  of  amplitude  is  not  useful  for 
discriminating  between  different  sets  of  similar  mechanisms. 
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V’avefom  Studies  of  the  Velocity  and  Q  Structure  of  the  Mantle 


TECHNICAL  SUMMARY 


In  recent  years,  the  techniques  of  time  domain  waveform  synthesis 
have  permitted  us  to  understand  much  more  about  seismic  sources,  about  the 
seismic  velocity  structure  of  the  earth  and  about  seismic  attenuation  than 
ever  before.  In  a  second,  separate  program  of  study,  which  is  described 
here,  an  effort  is  being  made  to  apply  these  techniques  to  solution  of  the 
detection,  discrimination  and  yield  estimation  problems.  A  variety  of 
investigations  have  either  been  completed  or  are  under  way  which  are 
designed  to  derive  as  much  information  about  the  earth  as  possible  from 
seismic  body  waves. 

We  present  the  ^^nal  results  of  a  study  of  the  structure  of  the  upper 
mantle  beneath  North  America  and  Europe.  Models  of  the  velocity  structure 
of  the  earth  were  developed  which  explain  the  major  features  of  P  wave 
propagation  at  ranges  between  IS*  and  30°.  Final  results  are  presented  on 
a  study  of  the  Gazli  earthquakes  in  Uzbekistan,  USSR.  They  occured  in  a 
region  which  is  geographically  convenient  for  sampling  wave  propagation 
within  the  Soviet  Union.  Consequently,  they  have  already  been  studied  by  a 
number  of  investigators  working  on  discrimination  problems  (Pearce,  1980; 
Pomeroy,  1980).  The  Gazli  area  had  been  relatively  aseismic  in  recent 
times,  but  an  analysis  of  the  historic  record  shows  that  it  is  very  pro¬ 
bably  a  seismic  zone  with  a  long  repeat  time.  A  waveform  analysis  of  the 
teleseismic  body  waves  is  also  presented.  Finally,  interim  reports  are 
presented  on  investigations  in  progress  of  a  major  three  dimensional 
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seismc  feature  of  geographic  importance  to  monitoring  sources  in  the 
Soviet  Union  and  of  the  frequency  dependence  of  Q  m  the  earth.  In  the 
former  study,  an  analysis  of  the  particle  notion  of  WWSSb  stations  located 
on  the  collision  boundary  between  the  Indian  and  Eurasian  continents  lias 
indicated  that  the  plate  boundary  has  a  significant  effect  on  long  period 
oody  wave  propagation.  Further  work  is  in  progress  to  explain  the  pheno¬ 
menon  in  terms  of  a  tectonic  model.  In  the  latter  study,  investigations 
are  being  undertaken  to  determine  how  significant  the  frequency  dependence 
of  Q  wi 1 1  be  to  waveform  analysis  techniques  and  also  to  determine  what  the 
frequency  dependence  of  Q  is. 

A .  A  Comparison  of  the  Upper  Mantle  Structure  Beneath  the  US  and  Europe 

The  characteristics  of  body  wave  propagation  to  epicentral  distances 
of  20°  or  more  depends  to  a  large  extent  on  tectonic  regime.  The  western 
US,  eastern  US  and  Europe  are  important  test  areas  for  understanding  this 
dependence  because  they  are  well  sampled  by  sources  and  stations  and 
because  tney  have  already  been  extensively  studied  by  many  researchers. 
Some  of  the  more  significant  investigations  include  the  travel  time 
studies  of  Hales  (19721  and  England  and  Worthington  (1477);  the  apparent 
velocity  studies  of  King  and  Calcagnile  11475),  England,  Kennett  and 
Worthington  (1973)  and  Johnson  (1907);  and  the  waveform  analysis  studies 
of  Wiggins  and  Helmberger  (1973),  Dey  Sarkar  and  Wiggins  (1975)  and  Given 
and  Helmberger  (1980).  Each  method  for  determining  mantle  structure  has 
provided  some  important  information,  but  generally  speaking,  travel  time 
studies  have  provided  the  roughest  picture  usually  showing  that  tectoni¬ 
cally  active  regions  are  slow  with  respect  to  stable  continent.  Apparent 
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velocity  studies  have  led  to  more  refined  models  and  wave  form  analysis 
studies  have  resolved  the  finest  detail.  The  new  work  presented  here  shows 
that  all  the  different  types  of  information  need  to  be  taken  into  account. 
Even  though  travel  time  data  provides  only  rough  constraints,  it  can  be 
important.  A  complete  report  on  the  upper  mantle  study  of  North  America 
and  Europe  is  presented  in  Appendix  I.  It  has  been  accepted  for  publi¬ 
cation  in  the  Journal  of  Geophysical  Research. 

In  a  parallel  study  to  the  one  reported  here  Burdick  and  Helmberger 
(1978)  completed  an  upper  mantle  study  of  the  western  US.  They  proceded  by 
first  modeling  several  seismic  sources  very  accurately  and  then  by 
modeling  the  structure  of  the  upper  mantle  using  records  from  the  known 
sources.  The  main  new  result  of  that  work  was  the  specific  identification 
in  long  period  body  waves  of  pulses  reflected  from  the  upper  mantle 
discontinuities.  Figure  1  shows  a  particularly  good  example  of  a  reflec¬ 
tion  (labeled  01)  from  the  670  km  discontinuity  in  the  mantle.  In  the 
course  of  their  study  they  developed  model  T7  which,  as  shown  in  the 
figure,  had  some  advantages  over  a  previous  model  HWA.  The  data  set  used 
to  arrive  at  model  T7  did  have  a  shortcoming,  however,  in  that  no  obser¬ 
vations  were  available  from  the  back  branch  of  the  670  km  discontinuity. 
The  back  branch  of  the  400  km  discontinuity  was  clearly  identified  in 
records  from  near  15°  and  its  forward  branch  in  records  from  near  23°.  The 
forward  branch  of  the  670  km  discontinuity  was  clearly  identified  in 
records  near  27°  (Figure  I),  but  because  appropriate  records  were  not 
available  from  near  19°,  model  T7  was  unconstrained  with  respect  to  the 
behavior  of  the  forward  branch  of  this  triplication.  Figure  2  shows 
several  records  from  the  new  data  set  at  ranges  of  19°  to  20°  which  do 
contain  a  small  feature  associated  with  the  missing  triplication  branch. 
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Long  period  records  from  the  ranges  26°  to  30c 
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Figure  1 


The  long  period  records  shown  in  the  figure  contain  two 
arrivals.  The  second  arrival  which  is  marked  as  D1  is 
a  reflection  from  the  670  km  discontinuity  (see  Figure 
3.4).  The  effective  source  function  is  shown  on  the  left 
for  each  station.  The  filtered  delta  function  response 
is  next  and,  the  synthetic  seismogram  is  on  the  right  in 
light  line.  The  observed  record  is  on  the  right  in  dark 
line.  Model  T7  appears  to  correctly  predict  the  behavior 
of  the  reflection,  but  HVA  does  not. 
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Figure  2  Figure  arrangement  is  the  same  as  figure  1.  A  further 
discussion  is  given  in  appendix  I. 
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Model  T9  is  Che  new  upper  mantle  model  introduced  in  Appendix  I  to  be 
consistent  with  the  new  data.  The  model  fits  the  arrival  at  19.3°  and 
i9.6°  and  as  discussed  in  the  appendix  can  also  fit  the  arrival  at  19.3°  it 
a  different  source  is  assumed.  Figure  3  shows  the  two  new  upper  mantle 
models  T9  (tectonic  continent)  and  S8  (stable  continent).  They  appear  to 
correctly  predict  the  behavior  of  long  period  P  wave  propagation  out  to 
ranges  of  30°.  As  has  long  been  suspected  from  surface  wave  studies,  the 
main  difference  becween  tectonic  and  stable  continent  appears  to  be  simply 
the  appearance  or  disappearance  of  a  low  velocity  zone. 

The  second  important  new  result  of  the  new  study  is  that  there  is  a 
baseline  difference  in  travel  time  between  models  for  Northeast  Eurasia 
(King  and  Calcagnile,  1976;  Given  and  Helmberger,  1980)  and  models  T9  and 
S8.  The  former  models  are  tied  to  ISC  origin  times  for  Russian  nuclear 
blasts,  and  the  latter  are  tied  to  NTS  origin  times.  Figure  U  compares  the 
model  cf  Given  and  Helmberger  (  1980)  K8  to  model  S3.  It  is  clear  that  most 
of  the  difference  in  travel  time  has  been  aosorbed  in  a  difference  in  depth 
to  Che  major  discontinuities.  Other  than  this  difference,  the  two  models 
are  very  similar.  The  ISC  times  for  the  Soviet  nuclear  tests  are  generally 
several  seconds  before  the  minute.  If  this  is  a  systematic  error,  then 
perhaps  the  depth  to  the  major  discontinuities  does  not  vary  as  signifi¬ 
cantly  as  would  be  indicated  by  Figure  4.  In  any  case,  it  would  now  appear 
that  in  the  models  of  Given  and  Helmberger  (1980),  in  T9  and  in  S8  we  have 
developed  very  useful  tools  for  understanding  the  properties  of  P  wave 
propagation  in  the  ranges  10°  to  30°. 
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B .  The  Focal  Mechanisms  of  the  Gazli  Earthquakes 


The  Gazli  sequence  of  earthquakes  in  Uzbekistan,  USSR  is  a  very 
important  one  because  of  its  unique  location.  No  comparable  sequence  has 
occurred  since  the  advent  of  modern  seismic  recording.  Pearce  (1980)  found 
the  sequence  to  be  very  useful  in  testing  his  technique  for  source  discrim¬ 
ination  using  array  data.  Pomeroy  (1980)  also  used  the  sequence  in  a  study 
of  Lg  propagation  in  the  Soviet  Union. 

The  Gazli  sequence  contained  two  major  shocks,  both  of  magnitude  7.1. 
Because  of  the  overall  importance  of  the  sequence,  we  have  undertaken  a 
waveform  analysis  study  of  these  two  events  to  determine  their  focal  mech¬ 
anisms.  We  have  also  studied  the  historical  seismicity  records  to  attempt 
to  determine  why  they  happened  where  they  did.  A  complete  report  is  given 
in  Appendix  II.  It  has  been  published  in  the  Bulletin  of  the  Seismological 
Society  of  America.  Figure  5  shows  the  locations  of  the  Gazli  earthquakes. 
Also  shown  as  a  shaded  area  is  an  arcuate  region  which  underwent  severe 
deformation  in  Hercynian  times.  The  Gazli  events  appear  to  have  occurred 
on  this  belt.  This  would  support  the  suggestion  of  Sykes  (1978)  that 
intraplate  earthquakes  tend  to  occur  on  such  belts.  Figures  6  and  7  show 
the  final  focal  mechanism  solutions  for  the  two  events  and  the  associated 
waveform  fits.  The  most  important  observation  is  that  the  two  events  had 
substantially  different  mechanisms.  It  is  often  assumed  that  the  large 
events  in  a  sequence  all  occur  on  the  same  fault  plane.  As  discussed  in 
Appendix  II,  Aptekman  et  al.  (1978)  have  suggested  that  this  may  not  be 
true  for  thrust  events.  The  evidence  from  Gazli  tends  to  bear  them  out. 
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C .  A  Major  Seismic  Feature  of  Geographic  Importance  to  Monitoring  Soviet 

Seismic  Events 


It  has  long  been  suspected  that  large  scale  3-dimensional  features 
such  as  subduction  zones  might  affect  our  ability  to  detect  and  character¬ 
ize  seismic  events.  Ordinarily,  important  seismic  installations  are  not 
situated  on  such  features  in  an  effort  to  avoid  the  problem.  In  spite  of 
these  efforts,  however,  LASA  and  NORSAR  were  inadvertently  located  at  less 
than  ideal  sites  (Aki,  1973;  Serteussen  et  al.,  1975).  burdick  and 
Langston  (1977)  and  Langston  (1977,  1978)  have  shown  that  many  WKSSN 
stations  are  strongly  affected  ay  3-dimensional  structure  at  the  receiver. 
A  survey  has  been  performed  of  the  site  characteristics  of  WVSSN  stations 
in  India  and  southern  Eurasia  to  identify  those  which  are  transparent.  In 
the  course  of  the  work,  it  was  found  that  stations  located  near  the 
collision  boundary  between  India  and  Eurasia  are  all  situated  on  strongly 
dipping  structure.  This  has  major  implications  for  the  monitoring  effort 
because  these  stations  and  the  new  digital  stations  which  are  replacing 
many  of  them  are  geographically  ideal  for  recording  Novaya  Zemlya  and 
Kazakh  events. 

A  survey  of  station  site  characteristics  must  generally  proceed  in 
two  parts.  First,  the  particle  motion  characteristics  of  the  stations  are 
examined  to  determined  if  three  dimensional  velocity  variations  are 
affecting  the  incoming  waves.  Then,  in  those  instances  where  the  velocity 
structure  appears  to  be  flat  lying,  plane  layered  models  for  the  crust  are 
constructed  using  P  to  S  converted  phases  in  teleseismic  body  waves.  The 
stations  considered  in  this  survey  are  shown  in  Figure  8.  The  shaded 
region  approximately  shows  the  extent  of  the  collision  boundary  between 
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India  and  Eurasia  (Saeber  et  al.,  1961).  The  stations  located  on  the 
boundary  are  QL'E,  K6L  ana  ML.  It  is  important  to  note  that  QUE  is  going  to 
be  upgraded  to  a  digital  WWSSL  station  and  K.EL  already  has  r.n  SRC. 

If  a  teleseisir.ic  (  A  >30°)  P  wave  impinges  on  a  flat  lying  crustal 
column  then  the  compress ional  motion  will  not  couple  to  the  tangential 
notion.  This  implies  that  the  north  component  of  motion  must  be  just  a 
constant  times  Che  east  component. 

U(t)  =  c  E(T) 


where  c  =  cot  (3ack  Azimuth). 

At  long  tines  (t  *  30  sec),  energy  scattered  far  from  the  station  will 
cause  this  relationship  to  break  down.  At  short  times,  however,  serious 
breakdown  of  the  relationship  can  only  be  caused  by  dipping  structure  at 
the  station.  The  quality  of  a  site  can  be  determined  by  simply  comparing 
the  two  horizontal  components  of  motion  for  teleseismic  P  waves  to  each 
other . 

Figure  9  shows  some  typical  horizontal  P  wave  components  for  stations 
in  the  survey.  The  column  on  the  left  shows  some  of  the  best  stations;  the 
column  in  the  center,  some  intermediate  stations  and  the  column  on  the 
rignt  the  very  poor  ones  which  are  those  on  the  collision  boundary.  The 
tangential  arrivals  at  the  poor  stations  are  very  impulsive  indicating  a 
sharp  velocity  increase  on  a  strongly  dipping  interface.  It  is  not 
uncommon  to  find  that  the  particle  motion  at  any  station  is  anomalous  from 
certain  azimuths.  Figures  10  and  11  show  additional  comparisons  of  hori¬ 
zontal  components  of  motion  at  KBl  and  QUE  (NIL  has  a  relatively  short 
recording  historv).  The  particle  motion  is  poor  for  all  azimuths  of 
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approach.  The  P  waves  coning  from  almost  due  north  (Back  Azimuth  =  353.6 
at  QU2  and  Back  Azimuth  =  353.7  at  KBL)  are  particularly  interesting.  They 
originated  from  a  Soviet  nuclear  test  at  h'ovaya  Zemlya.  Because  the 
station  is  nearly  naturally  rotated  the  east  component  is  almost  identical 
to  the  tangential  or  "SH"  component.  JJote  how  the  P  wave  bleeds  onto  the 
tangential  component  with  tine.  This  phenomenon  is  currently  being 
investigated  further,  but  from  just  this  survey  we  can  conclude  that  the 
continental  collision  boundary  is  a  major  three  dimensional  velocity 
structure  which  affects  seismic  recording  in  an  adverse  way. 

The  approach  used  in  the  second  part  of  the  survey  to  model  the 
crustal  structure  at  the  good  and  intermediate  quality  sites  is  the  one 
developed  by  Burdick  and  Langston  (1977).  The  P  to  S  converted  phases  in 
teleseisms  are  specifically  identified  and  modeled  using  synthetic  seismo¬ 
grams.  Figure  12  shows  theoretical  seismograms  for  a  plane  wave  incident 
from  below  on  a  layer  over  a  half-space.  The  vertical  component  of  motion 
is  shown  on  the  left  and  the  horizontal  on  the  right.  The  first  order 
crustal  reverberations  are  named  and  numbered  using  the  notation  of  Bath 
and  Steffanson  (1966).  The  S  phases.  A,  5  and  6,  cause  clear  differences 
between  the  two  components  of  motion.  The  largest  arrival  of  the  three  is 
tne  first  one,  the  Ps  phase.  We  have  been  able  to  observe  this  phase  in 
data  from  all  of  the  good  and  intermediate  quality  stations  in  Figure  9. 
3y  modeling  it,  we  have  been  able  to  constrain  the  crustal  thickness  in 
each  case.  Figure  13  compares  data  to  synthetics  for  station  SHI.  The  Ps 
arrival  is  clear  in  the  radial  components  and  is  well  matched  by  the 
synthetics.  The  crustal  thickness  is  33  km,  very  typical  for  continental 
crust.  (We  assume  throughout  that  =  8.1  km/ sec  in  the  lid.)  Similar 
results  for  stations  MSH,  ROD,  POO  and  NDI  are  shown  in  Figures  1A-17. 
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values  i  ;r  crustai  thickness  ranges  from  30  to  30  kn.  There  does  not 
appear  to  be  any  systematic  difference  between  the  Indian  and  Eurasian 
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D .  Relative  Attenuation  of  P  and  S  Waves  and  the  Frequency  Dependence  of 


a 


The  recent  attempts  to  correlate  the  results  of  body  wave  attenuation 
studies  over  a  broad  range  of  frequencies  have  produced  evidence  that  Q 
must  vary  with  frequency  (Sipkin  and  Jordan,  1S79).  This  raises  two 
important  questions.  First,  how  does  this  affect  the  results  of  studies 
which  had  assumed  a  constant  Q  and  second  what  is  the  frequency  dependence 
of  Q?  Both  of  these  questions  are  currently  being  investigated. 

Some  of  the  studies  which  have  assumed  that  Q  is  independent  of 
frequency  have  been  the  long  period  body  wave  modeling  studies  of 
Helmberger  (1S74),  Burdick  and  Mellman  (1976)  and  Langston  (1976).  It  is 
important  to  establish  whether  or  not  these  results  would  change  if  Q  were 
assumed  to  depend  on  frequency.  The  procedure  used  in  these  investigations 
to  correct  for  Q  was  to  include  a  Futterman  (1962)  constant  Q  operator  in 
the  calculation  of  the  synthetic  seismograms.  Sipkin  and  Jordan  (1979) 
have  suggested  that  an  appropriate  modification  would  be  to  use  a  standard 
linear  solid  (S.L.S.)  Q  operator  (Minster,  1978).  The  S.L.S.  operator  has 
three  adjustable  parameters.  One  is  just  the  minimum  value  of  Q  which  is 
the  equivalent  of  t*  for  a  Futterman  operator.  The  other  parameters 
control  the  increase  of  Q  at  high  and  low  frequencies.  The  low  frequency 
behavior  of  Q  is  not  important  to  body  wave  calculations,  but  the  high 
frequency  behavior  is.  The  controlling  parameter  is  t  .  Sipkin  and  Jordan 
(1979)  suggest  that  it  should  be  between  .1  and  1.0  sec,  but  as  we  show  in 
the  following  section  it  is  probably  at  the  low  end  of  this  range. 


45 


Figure  IS  is  adapted  from  the  work  of  Sipkin  and  Jordan  and  summarizes 

their  arguments.  The  hatched  areas  show  esimates  of  the  average  shear  Q  of 

the  mantle  (Qg)  as  determined  from  normal  modes,  ScS^  experiments  and 

spectral  decay  measurements.  It  would  seem  reasonable  on  the  basis  of  this 

evidence  to  adont  a  value  of  t  =  .5.  However,  careful  consideration  of 

m 

the  results  of  measurements  of  the  relative  attenuation  of  P  and  S  waves 
along  raypaths  where  the  source  excitation  is  nearly  constant  show  that 
this  is  not  the  case. 

We  begin  by  noting  that  the  high  frequency  spectral  decay  measure¬ 
ments  of  Sipkin  and  Jordan  are  much  less  reliable  then  their  low  frequency 
ScS^  measurements.  The  spectral  decay  measurements  are  made  by  solving  for 
the  attenuation  spectrum  A(u)  from  an  observation  S(uj)  using 


S(w)  =  A (w)  R(u> ) 


A(uj).  the  spectrum  of  the  source  is  also  unknown.  The  large  vertical  error 
bars  on  the  spectral  decay  measurements  in  the  figure  were  derived  by 
testing  a  range  of  source  models.  The  dark  horizontal  bars  are  based  on 
the  source  parameter  measurements  of  Wyss  and  Molnar  (197  2).  Wyss  and 
Molnar  also  solved  an  equation  like  (2),  but  they  assumed  a  form  for  A(;u) 
and  solved  for  R(oj).  Their  Q  correction  was  based  on  the  results  of  Julian 
and  Anderson  (  1968)  who  used  model  MM8.  The  Qg  for  this  model  is  580  which 
is  very  close  to  the  value  recovered  by  Sipkin  and  Jordan.  That  is  to  say, 
the  source  spectrum  and  attenuation  spectrum  determinations  are  sel  f- 
consistant  but  totally  non-unique.  This  does  not  invalidate  Sipkin  and 
Jordan's  arguments  for  a  lower  bound  for  Q  because  seismic  sources  must 
have  corner  frequencies.  It  does  show,  however,  how  important  it  is  to 
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make  measurements  of  attenuation  in  situations  where  the  source  excitation 
is  controlled. 

One  way  to  measure  attenuation  so  that  the  source  excitation  is 
controlled  is  to  measure  the  relative  attenuation  of  ScS  with  respect  to 
SoP.  In  doing  so,  Kanamori  (  1967)  measured  Qg  for  the  mantle  to  be  between 
150  and  260  in  the  frequency  range  of  .1  to  .7  Hz.  Sacks  (1980)  repeated 
the  experiment  and  obtained  a  value  of  250  between  .3  and  1.5  H2 .  Burdick 
(1978)  compared  sP  to  sS  phases  and  obtained  a  value  of  tg*  =  5.2  *_  .7  sec 
at  a  range  of  A  =  40°.  The  5L8  Q  model  (Anderson  and  Hart,  1979)  predicts 
tg*  =  4.6  sec  at  this  distance  and  has  an  overall  of  235.  The  Burdick 
(1978)  measurement  therefore  scales  to  210  +  25  for  WWSSN  short  period 
records  low  pass  filtered  at  1  Hz. 

The  relative  S-P  attenuation  measurements  are  represented  in 
Figure  19  as  a  stipled  region  of  Q  bounded  between  150  and  260  over  a 
frequency  range  of  .2  to  1.5  Hz.  The  measurements  clearly  suggest  that  the 
rapid  increase  in  Q  with  frequency  does  not  begin  to  occur  until  after 
l  Hz.  The  cutoff  time  t  might  be  appropriately  chosen  as  .i  to  .2  sec 
which  reconciles  as  closely  as  possible  the  observations  of  Sacks  (1930) 
with  the  lower  bound  of  Sipkin  and  Jordan  (1979).  This  would  imply  that  Q 

□ 

rises  to  1000  by  5  Hz.  An  appropriate  constant  Qg  approximation  for  all 
the  data  in  the  WWSSN  band  (.01  to  2  Hz)  would  be  a  value  of  250.  The 
important  point  here  is  that  measurements  of  the  relative  attenuation  of  P 
and  S  waves  in  cases  where  the  source  excitation  is  controlled  give  a  very 
self-consistent  result.  They  appear  to  show  that  %  does  noc  rise  signifi¬ 
cantly  above  300  at  frequencies  less  than  1.5  Hz. 

To  determine  the  importance  of  the  frequency  dependence  of  Q  to  long 
period  waveform  studies,  synthetic  seismograms  for  the  Borrego  Mountain 
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earthquake  (Surdick  and  Mellman,  1976)  were  computed  using  a  Futtermar.  and 

an  S.L.S.  operator  (t  =  .2  sec).  The  results  are  compared  in  Figure  19. 

m 

The  observed  P  wave  from  COL  and  the  observed  SH  from  OGD  are  shown  at  the 
top.  Underneath  them,  superimposed  on  each  other,  are  synthetics  computed 
with  the  two  different  operators.  The  differences  are  almost  invisible. 
From  this  and  from  other  tests,  we  have  concluded  that  the  fact  that  Q  does 
increase  with  frequency  will  have  no  effect  on  the  results  of  the  long 
period  waveform  modeling  studies. 
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CONCLUSIONS 


The  studies  of  the  Gazli  earthquakes  and  of  the  upper  mantle  show  that 
we  are  becoming  very  proficient  at  predicting  many  of  the  features  of  long 
period  body  waves  using  simple  point  source  approximations  and  plane 
layered  earth  structures.  These  approaches  will  undoubtedly  continue  to 
be  of  value  to  us.  The  results  of  the  station  survey,  on  the  other  hand, 
show  that  three  dimensional  velocity  structures  do  have  a  significant 
effect  on  long  period  body  waves  in  some  instances.  The  effect  will 
presumably  be  even  larger  on  short  period  body  waves.  Likewise,  the 
frequency  dependence  of  Q  will  be  important  to  the  analysis  of  short  period 
body  waves.  Since  many  small  events  are  recorded  only  at  short  periods  it 
will  be  important  to  continue  our  work  in  these  areas. 
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Appendix  I 


A  COMPARISON  OF  THE  UPPER  MANTLE  STRUCTURE 


BENEATH  NORTH  AMERICA  AND  EUROPE1 

L.  J.  Burdick 

Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  New  York  10964 

ABSTRACT  • 


The  techniques  of  modeling  upper  mantle  structure  by  matching  long 
period  waveforms  with  synthetic  seismograms  have  been  applied  to  obser¬ 
vations  from  the  tectonically  stable  part  of  North  America  and  from  Europe. 
The  consistent  differences  which  can  be  resolved  by  the  long  period  data 
between  Europe  and  North  America  can  be  interpreted  in  terms  of  variations 
in  the  crust,  lid  and  low  velocity  zone.  At  epicentral  ranges  less  than 
15°,  the  effects  of  shallow  lateral  variations  are  strong  and  body  wave 
propagation  is  regionally  dependent.  Between  ranges  of  15°  and  20°, 
regional  effects  are  still  observed,  but  they  can  be  explained  in  terms  of 
variations  above  250  km.  Beyond  20°,  wave  propagation  appears  to  be  stable 
and  independent  of  region.  Most  of  the  observed  long  period  P  waves  from 
20°  to  30°  are  consistent  with  a  single  model.  This  indicates  that  the 
relative  depths  and  sizes  of  the  major  discontinuities  do  not  vary  sub¬ 
stantially.  A  comparison  of  the  upper  mantle  models  of  this  study  with 
those  of  other  studies  indicates  that  the  shape  of  the  P  velocity  profile 
is  fairly  uniform  through  the  transition  region  (300-700  km)  though  there 
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are  differences  in  the  absolute  depth  to  the  discontinuities.  These 
differences  may  be  the  result  of  systematic  errors  in  travel  time  data. 


INTRODUCTION 


Much  recent  interest  in  the  problem  of  lateral  variations  in  mantle 
structure  has  been  caused  by  the  suggestion  of  Jordan  (1975)  that  the 
structural  roots  of  continents  extend  down  to  400  km  or  more.  '  The  primary 
support  for  this  suggestion  comes  from  studies  of  ScS  travel  times  and 
surface  wave  dispersion  [Sipkin  and  Jordan,  1975,  1976;  Jordan,  19791. 
Opposing  evidence  showing  that  continents  are  much  shallower  than  400  km 
has  been  presented  by  Okal  and  Anderson  [19751  and  Okal  [1977,  19781.  The 
investigations  of  higher  mode  Rayleigh  waves  by  Cara  [19781  and  of  P  wave 
travel  times  and  apparent  velocities  by  England  et  al.  [19781  also  indicate 
that  major  structural  differences  between  continents  and  oceans  are  con¬ 
fined  to  depths  above  250  km.  The  question  is  still  very  much  under 
debate.  An  issue  that  is  just  as  important  in  resolving  the  depth  of 
continental  roots  is  the  depth  of  lateral  variations  between  and  within 
continents  themselves  [Jordan,  19791.  Both  Sipkin  and  Jordan  [19761  and 
England  £t  al_.  [19771  have  found  evidence  that  structural  differences 
within  continents  extend  well  below  250  km.  At  the  current  time,  there 
appears  to  be  at  least  as  much  evidence  to  support  the  idea  of  deep  lateral 
variations  between  continents  as  there  is  to  support  variations  between 
continents  and  oceans. 

An  unambiguous  picture  of  the  extent  and  magnitude  of  lateral  vari¬ 
ations  in  seismic  velocity  throughout  the  upper  mantle  will  most  certainly 
help  in  mapping  the  course  of  its  geochemical  and  tectonic  evolution. 
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However,  such  a  picture  is  hard  to  derive  from  comparisons  of  the  results 
of  different  types  of  studies  in  different  regions  of  the  world.  It  13 
difficult  to  sort  out  true  variations  in  structure  from  apparent  vari¬ 
ations  caused  by  differences  in  methodology.  Here  we  present  the  results 
of  an  effort  to  compare  the  P  velocity  structures  of  North  America  and 
Europe  in  the  depth  range  of  250  to  1000  km  using  the  single  method  of 
modeling  waveforms  with  synthetic  seismograms. 

The  problem  of  resolving  large  scale  lateral  variations  in  the  250  km 
to  1000  km  depth  range  is  severely  complicated  by  the  fact  that  there  are 
substantial  variations  above  250  km  in  the  crust,  lid  and  low  velocity 
zone.  Waveform  data  usually  come  from  sources  in  tectonically  active 
regions  and  are  recorded  at  stations  in  both  active  and  stable  regions. 
The  process  of  distinguishing  deep  variations  in  structure  from  shallow  is 
a  very  difficult  one.  In  this  investigation,  we  have  attempted  to  find  the 
distance  range  over  which  clear  evidence  for  lateral  variation  in  struc¬ 
ture  persists.  We  have  attempted  to  find  whether  observed  variations  in 
wave  propagation  can  be  explained  in  terms  of  shallow  structural  differ¬ 
ences  alone  or  whether  deeper  differences  are  indicated.  Finally  we  have 
attempted  to  define  a  set  of  models  which  can  explain  the  different  wave 
propagation  characteristics  of  P  waves  within  the  regions  studied. 


SYNTHETIC  SEISMOGRAM  METHODS 


The  techniques  for  modeling  upper  mantle  velocity  structure  using 
synthetic  seismograms  have  been  described  and  used  by  Wiggins  and 
Helmberger  119731,  Dey  Sarkar  and  Wiggins  (19761  and  by  Burdick  and 
Helmberger  (19781.  The  procedure  consists  of  two  parts.  First  accurate 
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source  models  for  moderate  sized  earthquakes  or  explosions  are  determined 
using  data  from  ranges  where  earth  structure  does  not  significantly  affect 
the  wave  shape  (30s  to  90s).  Then  synthetics  are  computed  to  match  wave¬ 
form  observations  from  ranges  where  upper  mantle  structure  has  a  strong 
effect  (15s  to  30s).  A  starting  model  is  chosen  for  the  upper  mantle 
structure,  and  its  response  for  a  given  range  is  computed  using  the 
Cagniard-de  Hoop  technique  (Wiggins  and  Helmberger,  1974].  The  response 
is  convolved  with  the  source,  instrument  and  attenuation  operators  to 
produce  the  synthetic.  If  the  final  synthetic  does  not  match  the  data,  the 
starting  model  is  perturbed  by  trial  and  error  until  a  satisfactory  fit  is 
attained. 

A  more  detailed  discussion  of  the  methods  for  computing  synthetics  is 
unnecessary  in  light  of  the  previous  contributions.  However,  some  recent 
studies  have  raised  some  issues  that  need  to  be  mentioned.  One  possible 
difficulty  has  been  brought  to  light  by  the  investigations  of  Burdick  and 
Orcutt  [1979]  who  compared  the  results  of  the  Cagniard-de  Hoop  method  to 
those  of  the  reflectivity  method  for  upper  mantle  structures.  They  showed 
that  the  primary  ray  approximation  breaks  down  in  the  Cagniard-de  Hoop 
method  for  thin  transition  zones.  The  primary  ray  sum  is  adequate  only  for 
models  with  moderate  gradients  and  first  order  discontinuities.  The 
extended  ray  sum  which  is  necessary  for  thin  transition  zones  has  not  yet 
been  established  and  the  general  behavior  of  the  seismic  wave  field  in  the 
presence  of  high  gradients  is  not  well  understood.  Both  of  these  questions 
must  be  answered  by  future  research,  so  only  appropriately  simple  models 
will  be  discussed  here. 

A  second  point  which  has  received  some  attention  recently  involves 
the  accuracy  of  the  correction  for  attenuation.  In  most  cases  a  Futterman 


[1962]  operator  has  been  used.  New  work  by  Sipkin  and  Jordan  [1979]  among 
others  has  shown  that  Q  may  begin  to  increase  with  frequency  at  some  point 
within  the  seismic .band.  These  effects  are  most  important  for  short  period 
body  waves,  however,  and  here  we  will  be  concerned  with  long  period  data. 
Both  Kennett  [1975]  and  Burdick  and  Helmberger  [1978]  have  considered  the 
importance  of  the  depth  dependence  of  Q  on  synthetic  calculations.  The 
latter  showed  that  the  current  method  is  adequate  if  the  Q  distribution  in 
the  earth  falls  within  a  rather  broad  range.  However,  they  investigated  a 
specific  velocity  model  with  a  relatively  thin  lithosphere.  Rays  which 
traveled  in  this  high  Q  region  were  too  slow  to  be  of  significance  beyond 
10°.  In  a  model  with  a  thicker,  faster  lithosphere,  these  unusually  high  Q 
paths  become  much  more  important.  Such  regional  phases  as  PL  and  Pr  may 
contaminate  the  waveform  out  to  much  greater  ranges.  In  this  study,  the 
data  set  has  been  limited  to  observations  from  beyond  14*  to  minimize  these 
possible  problems. 


THE  REGIONS  OF  STUDY 


The  distribution  of  WWSSN  stations  and  active  seismic  zones  in  the 
world  make  the  U.S.  and' Europe  the  two  regions  that  are  most  suitable  for 
modeling  mantle  structure  with  the  synthetic  seismogram  approach.  The 
central  and  eastern  U.S.  can  be  studied  using  sources  in  the  intermountain 
seismic  belt  on  one  side  and  in  the  Caribbean  on  the  other.  Central  and 
southern  Europe  can  be  studied  using  sources  in  the  North  Aegean  region 
with  stations  in  Scandinavia  or  western  Europe.  Much  is  already  known 
about  the  basic  character  of  the  mantle  in  both  regions  because  careful 


studies  of  the  travel  times  have  been  made  for  each. 
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Hales  [1972]  has  summarized  the  work  on  travel  times  in  the  U.S. 
Observations  from  the  GNOME  explosion  in  eastern  New  Mexico  and  the  Early 
Rise  explosions  in  Lake  Superior  show  that  the  stable  part  of  the  North 
American  continent  is  very  fast  out  to  20°.  The  tectonically  active 
western  U.S.  is  as  much  as  6  seconds  slower  in  this  range  for  P.  Figure  1 
shows  the  locations  of  the  seismic  events  used  in  this  study  with  respect 
to  the  stations  in  the  U.S.  and  Canadian  networks.  The  first  event,  the 
3/28/75  Pocatello,  Idaho  earthquake,  occurred  on  the  eastern  edge  of  the 
intermountain  seismic  belt  [Smith  and  Sbar,  19741.  It  had  a  simple 
impulsive  source  which  produced  a  large  .  quantity  of  teleseismic  upper 
mantle  data.  Many  of  the  raypaths  from  the  source  to  the  eastern  and 
northeastern  stations  remain  exclusively  in  the  tectonically  stable  part 
of  the  continent.  It  should  be  noted  that  a  few  of  the  paths  to  stations  in 
the  south  or  northwest  (DAL,  JCT,  FSJ)  do  sample  the  tectonic  part  of  the 
continent  modeled  by  Burdick  and  Helmberger  [1978]  (Model  T7).  We  will 
show  that  along  these  paths  the  mantle  response  is  somewhat  different  than 
it  is  in  the  east.  The  second  event,  the  3/24/78  Bermuda  earthquake,  is 
also  well  located  for  sampling  the  seismic  structure  of  the  eastern  part  of 
the  continent.  In  this  case,  of  course,  the  beginning  of  the  raypath  is 
always  in  an  oceanic  regime.  For  those  paths  actually  used  in  the  study, 
the  ray  bottoming  points  are  primarily  under  the  continent  or  the  contin¬ 
ental  shelf.  The  location  of  the  4/9/68  Borrego  Mountain  earthquake  in 
California  is  also  shown  in  Figure  1  since  data  from  it  will  also  be 
considered. 

Figure  2  shows  the  sources  and  stations  used  in  studying  the  mantle 
structure  of  Europe.  The  travel  times  for  this  same  combination  of  source 
regions  and  stations  were  studied  by  England  and  Worthington  [1977].  They 


found  that  the  times  were  generally  much  slower  out  to  ranges  of  20®  than 
they  are  in  central  and  eastern  North  America.  In  this  sense  the  region 
should  perhaps  be  classified  as  tectonically  active  continent,  but  some  of 
the  WWSSN  stations  to  the  north  are  situated  in  areas  which  have  been 
stable  since  Hercynian  times.  (In  classifying  areas  as  tectonically 
active  or  stable  we  will  generally  adopt  the  map  given  by  Jordan  [ 19751) • 
So,  as  in  the  North  American  study,  some  of  the  raypaths  traverse  mixed 
tectonic  regimes.  The  two  North  Aegean  earthquakes  selected  for  use  are 
moderate  sized  normal  events  which  are  very  comparable  to  the  Idaho  event 
in  North  America.  They  produced  simple  clear  P  wave  records  at  all  ranges 
from  15®  to  90°. 


THE  SOURCE  MODELS 


The  earthquake  sources  are  modeled  as  simple  point  shear  dislocations 
buried  in  a  layered  half  space.  The  far  field  time  function  is  assumed  to 
be  a  triangular  pulse  specified  only  by  a  rise  and  fall-off  time  which  is 
independent  of  azimuth  and  take-off  angle  [Langston  and  Helmberger,  19751. 
The  final  source  models  based  on  fits  to  the  P  and  SH  waveforms 
from  ■  30®  to  90®  are  given  in  Table  1. 

The  Pocatello,  Idaho  earthquake  source  was  modeled  independently  by 
both  Burdick  [1978]  and  Bache  e£  al^.  [  19801  .  The  first  motion  data  alone 
constrained  the  fault  plane  solution  and  as  shown  in  Figure  3,  the  source 
model  fits  the  WWSSN  records  from  the  greater  distances.  An  accurate 
source  model  for  the  Bermuda  earthquake  was  first  presented  by  Stewart 
[1979].  (The  model  used  here  was  modified  slightly  to  fit  additional  first 
motion  data  from  South  America.)  It  was  a  thrust  earthquake  which  occurred 
under  about  5  km  of  water.  The  fits  of  the  P  waveforms  are  shown  in 
Figure  4.  The  closely  matched  details  later  in  the  waveform  are  caused  by 
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reverberations  in  tne  water  layer.  The  two  North  Aegean  events  were 
previously  studied  by  McKenzie  [1972]  using  first  motion  data.  His  pub¬ 
lished  source  mechanisms  were  not  accurate  enough  to  fit  the  waveform 
observations,  so  some  modifications  were  necessary.  The  new  solutions 
still  represent  dipping  normal  faults,  but  the  azimuths  of  the  planes  and 
the  amount  of  strike-slip  motion  are  different.  Both  P  and  SH  waveforms 
were  used  to  constrain  the  solutions.  The  final  waveform  fits  are  shown  in 
Figures  5  and  6. 

The  quality  of  fit  of  synthetics  to  data  in  Figures  3  through  6  should 
give  some  measure  of  the  accuracy  with  which  the  effects  of  the  source  and 
the  instrument  on  the  waveforms  can  be  predicted.  Though  the  fit  is 
generally  good,  it  is  not  exact.  The  degree  of  misfit  of  synthetics  at 
upper  mantle  ranges  must  be  compared  to  the  misfit  from  these  greater 
ranges  in  deciding  whether  it  results  from  failure  of  the  upper  mantle 
model  or  failure  of  the  source  model. 

MODELS  FOR  SHALLOW  VARIATION  IN  STRUCTURE 


The  problem  of  resolving  deep  lateral  variations  in  structure  in  the 
presence  of  strong  shallow  variations  is  a  difficult  one.  It  is  made  even 
more  difficult  by  the  fact  that  the  synthetic  seismogram  methods  to  be  used 
are  appropriate  for  purely  radial  structures.  The  simple  approach  we  have 
taken  is  to  define  a  model  that  we  take  to  be  an  average  for  tectonically 
active  regions  and  a  model  that  is  an  average  for  stable  regions.  The 
differences  in  the  models  above  250  km  are  designed  to  satisfy  observed 
differences  in  travel  time  and  vertical  delay  and  to  be  qualitatively 
consistent  with  surface  wave  regionalization  results.  From  these  two 
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models,  we  have  attempted  Co  establish  the  magnitude  of  the  effect  of  the 
shallow  lateral  variations  on  long  period  P  wave  propagation,  and  the 
distance  range  to  which  it  extends.  We  have  then  compared  the  results  with 
the  observed  properties  of  wave  propagation  in  North  America  and. Europe  to 
see  if  variations  occur  which  fall  outside  of  these  bounds. 

The  models  tested  are  shown  in  Figure  7  and  listed  in  Table  2. 
Tectonic  continent  model,  T9,  is  a  variant  on  the  western  U.S.  model  T7  of 
Burdick  and  Helmberger  119781.  (We  adopt  the  convention  of  showing  layers 
as  homogeneous  in  both  the  radial  and  flat-lying  case.)  The  only  differ¬ 
ence  is  that  the  lower  discontinuity  has  been  raised  slightly  for  reasons 
to  be  discussed  later.  The  model  has  a  well  developed  low  velocity  zone 
which  is  appropriate  for  a  tectonically  active  continent  [Johnson,  1967; 
Helmberger,  19731.  Stable  continent  model  S8  has  a  thicker  crust  and  no 
low  velocity  zone.  Below  250  km,  models  S8  and  T9  are  identical. 

Figure  8  compares  the  travel  time  curves  for  the  two  models.  The 
predominant  effect  of  the  low  velocity  zone  in  model  T9  is  to  slow  it  down 
with  respect  to  S8.  The  crosses  in  the  figure  are  travel  times  along 
stable  continent  paths  in  North  America.  They  include  the  Early  Rise  times 
along  the  Texas  and  Arkansas  profiles,  times  from  Early  Rise  to  some  other 
selected  stations,  times  from  Salmon  and  also  times  from  the  GNOME  explo¬ 
sion  along  the  eastern  profile  [Romney  al_ ,  ,  1962;  Jordan  et  al^.  ,  1966; 
Iyer  et_  £l.,  1969].  Green  and  Hales  (1968)  used  thi9  data  in  constructing 
their  stable  continent  model  ER1.  The  circles  are  travel  times  from 
tectonic  continent  paths  and  mixed  regime  paths.  They  include  the  data  set 
of  NTS  travel  times  considered  by  Burdick  and  Helmberger  11978]  and  the 
average  travel  times  for  Europe  given  by  England  and  Worthington  [1977). 
The  two  data  sets  clearly  define  the  travel  time  differences  between  stable 
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and  active  continent,  and  models  S8  and  T9  are  shown  to  be  representative 
of  the  two  types  of  behavior.  The  difference  in  vertical  delay  time  is 
1  second  between  the  two  models.  This  corresponds  with  the  measured 
average  residual  difference  between  the  eastern  and  western  U.S.  of  Carder 
et  al.  (19661.  The  qualitative  differences  between  the  models  are  consis¬ 
tent  with  the  well  known  results  of  surface  wave  regionalization  work 
[Knopoff,  19721.  The  differences  in  the  theoretical  responses  of  the 
models  should  provide  an  estimate  of  the  average  effects  of  shallow  differ¬ 
ences  in  structure  between  stable  and  tectonic  continent. 

Burdick  and  Helmberger  (19781  observed  that  at  distances  less  than 
15°  shallow  lateral  variations  in  structure  have  a  large  effect.  The 
characteristics  of  the  propagation  of  long  period  body  waves  are  so 
regionally  dependent  that  modeling  waveforms  with  an  average  earth  model 
is  not  productive.  The  same  is  true  at  short  ranges  in  the  regions 
considered  here.  The  effects  of  shallow  lateral  variations  are  less  severe 
beyond  15B,  but  they  exist  nonetheless.  Figure  9  compares  the  delta 
function  responses  for  S8  and  T9  at  ranges  between  15*  and  19*.  Even 
though  the  direct  ray  is  turning  below  250  km  and  the  large  reflection 
come3  from  400  km,  the  variations  above  250  km  have  a  significant  effect. 
The  relative  time  and  amplitude  of  the  two  arrivals  are  measurably  differ¬ 
ent.  We  can  expect,  therefore,  that  variations  in  lid  and  low  velocity 
zone  structure  will  continue  to  have  an  effect  in  these  ranges.  The 
magnitude  of  the  effect  will  typically  be  of  the  order  of  the  differences 
in  Figure  9,  but  because  S8  and  T9  are  average  rather  than  extremal  models 
it  could  be  even  larger.  Given  and  Helmberger  (1980)  show  that  if  the  lid 
velocities  are  slightly  higher  than  in  S8,  the  first  arrival  will  be  in  the 
lid  shadow  zone  out  to  19*.  It  will  be  very  fast  but  very  small.  Alter- 
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natively,  some  ray  paths  may  pass  through  an  even  more  substantial  low 
velocity  zone  than  the  one  in  T9,  so  that  the  two  arrivals  will  be  moved 
still  closer  together.  Thus,  if  we  find  in  the  data  variations  in  the 
relative  amplitude  and  time  of  the  two  arrivals  at  ranges  between  15s  and 
19°,  we  can  attribute  the  effect  to  lateral  varitions  above  250  km.  The 
interpretation  will  not  be  unique,  however,  and  we  cannot  preclude  vari¬ 
ations  from  250  to  400  km. 

By  19°  in  Figure  9,  the  responses  of  the  two  models  have  become 
similar.  The  lateral  variations  above  250  km  still  do  have  a  residual 
effect  beyond  this  range  because  the  distance,  4,  for  a  given  ray  parameter 
will  be  shifted.  Figure  8  shows  that  the  triplications  of  model  T9  are 
just  those  of  S8  shifted  to  shorter  range  by  about  half  a  degree.  In  real 
data  we  can  expect  to  sometimes  find  that  the  triplications  are  shifted  to 
different  ranges  for  different  raypaths  simply  because  of  lateral  vari¬ 
ations  above  250  km. 


THE  UPPER  MANTLE  TRIPLICATIONS 


The  two  discontinuities  in  models  S8  and  T9  cause  two  major  tripli¬ 
cations.  The  fir3t  begins  at  aDout  14°  and  ends  at  23°,  and  the  second 
extends  roughly  from  17*  to  26®.  Past  experience  has  shown  that  given 
branches  of  the  triplications  are  best  observed  over  selected  distance 
ranges.  These  are  the  ranges  where  the  triplication  branch  is  about  5  to 
10  seconds  behind  the  first  arrival.  Then  the  secondary  arrival  is  far 
enough  behind  the  first  so  that  it  is  clear  of  it  but  not  so  far  back  chat 
it  is  lost  in  the  noise  of  the  coda.  In  the  following,  we  will  consider  the 
data  in  groups  by  distance  range  according  to  which  triplication  branches 
have  the  largest  effect. 
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Figure  10  compares  data  to  synthetics  for  ranges  where  the  first 
arrival  comes  from  above  the  400  km  discontinuity,  and  a  strong  second 
arrival  is  caused  by  the  forward  branch  of  the  first  triplication.  As 
Figure  9  shows,  there  are  substantial  differences  between  $8  and  T9  at 
these  ranges,  so  in  Figure  10  synthetics  are  shown  for  each.  The  large 
two-sided  arrival  in  the  delta  function  responses  is  the  400  km  dis¬ 
continuity  reflection  and  the  small  one-sided  third  arrival  is  the  narrow 
angle  reflection  from  the  650  km  discontinuity.  The  MNT  and  OTT  records 
from  Bermuda  show  a  clear  feature  associated  with  the  triplication  arrival 
as  indicated  by  the  small  arrows.  Model  S8,  without  the  low  velocity  zone, 
fits  this  feature  much  better  than  T9.  All  of  the  Bermuda  records  are  more 
consistent  with  S8  indicating  that  this  model  is  more  appropriate  for 
eastern  North  America.  The  data  from  the  Idaho  event  are  not  more  con¬ 
sistent  with  one  model  than  the  other,  but  as  shown  in  Figure  1  many  of  the 
raypaths  traverse  partially  tectonic  and  partially  stable  areas.  The 
records  from  the  European  events  are  more  consistent  with  the  S8  model,  but 
again  the  raypaths  involved  traverse  mixed  regions. 

Evidence  for  some  type  of  lateral  variation  in  structure  is  abundant 
in  data  from  less  than  19°.  A  good  example  in  Figure  10  can  be  found  in  a 
comparison  of  the  records  from  DAL  and  JCT  for  the  Idaho  event.  Both 
stations  are  at  the  3ame  distance  (15.5°)  and  are  at  very  nearly  the  same 
azimuth  from  the  source  (Figure  1),  but  the  waveforms  are  different.  It 
will  be  shown  that  in  this  instance,  as  in  several  others,  the  discrepancy 
can  be  explained  by  a  slight  shift  in  the  arrival  time  of  the  reflected 
energy.  The  records  can  be  modeled  with  the  upper  mantle  response  of  T9  or 
S8  from  a  shifted  range  which  is  less  than  1®  away  from  the  true  range. 
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Between  19°  and  22®  there  are  five  critically  reflected  ray  arrivals 
causing  a  quintuplicati on.  Figure  11  compares  the  observations  to  the 
synthetics  for  these  ranges  from  model  S8.  (Model  T9  gives  almost  ident¬ 
ical  results.  1  The  delta  function  responses  show  the  strong  interaction 
between  the  various  arrivals.  The  fit  between  the  data  and  synthetics 
indicates  that  the  model  responses  are  accurate  at  these  very  sensitive 
ranges.  The  arrows  show  a  small  secondary  arrival  which  appears  at  19.3®, 
19.6"  and  19.8®.  This  feature  has  been  modeled  as  the  forward  branch  of 
the  second  triplication.  If  this  interpretation  is  correct,  these  are 
among  the  first  observations  of  this  branch  which  have  been  modeled  in  long 
period  waveform  data.  To  match  the  timing  of  this  phase,  it  was  necessary 
to  move  the  lower  discontinuity  up  to  650  km.  The  velocity  gradient 
beneath  the  discontinuity  was  reduced  to  compensate  for  the  change.  This 
keeps  the  model  responses  of  T7  and  T9  the  same  at  greater  distances.  The 
only  significant  differences  in  the  T7  and  T9  responses  occur  from  19°  to 
22®  where  the  arrival  time  of  the  lower  discontinuity  reflection  is  differ¬ 
ent.  Burdick  and  Helmberger  1  1 978 ]  did  not  have  any  observations  of  this 
reflection  so  T7  was  unconstrained  with  regard  to  this  triplication 
branch.  The  match  of  the  timing  and  amplitude  of  the  small  arrival  is  good 
at  19.3®  and  19.6®,  but  the  fit  at  MAI.  at  19  8®  is  relatively  poor. 
Nonetheless,  the  arrival  is  clear  in  the  data  and  some  possible  reasons  for 
the  misfit  are  discussed  in  the  following. 

The  data  from  ranges  between  15®  and  20®  contained  many  indications  of 
lateral  variation  in  mantle  response.  Figure  i2  shove  some  specific 
examples  in  which  the  fit  to  the  waveforms  can  be  much  improved  by  using 
the  theoretical  responses  from  a  shifted  range.  The  fir9t  and  second 
arrivals  vary  in  relative  timing  and  amplitude  by  about  a  second,  but  there 


is  no  evidence  for  anomalously  large  or  anomalously  small  second  arrivals 
and  no  evidence  for  additional  arrivals.  Also  shown  in  Figure  12  is  the 
record  from  19.8®  which  contains  an  arrival  from  the  lower  discontinuity. 
The  duration  of  the  source  time  function  in  the  synthetic  has  been  shor¬ 
tened  which  improves  the  fit  substantially.  This  indicates  that  source 
directivity  may  have  been  an  important  factor.  Attempts  to  improve  the  fit 
further  by  perturbing  the  upper  mantle  structure  were  unsuccessful. 
However,  the  amplitude  of  the  arrival  appears  to  be  nearly  correct,  and  at 
this  range  the  interference  between  arrivals  is  very  delicate.  A  3mall 
change  in  the  triplication  arrival  times  might  account  for  the  effect,  but 
a  major  change  in  the  size  of  the  discontinuity  is  not  suggested  by  the 
data.  Further  complications  in  the  source  or  structure  at  the  station  may 
also  have  had  an  effect. 

Figure  13  shows  data  and  synthetics  for  ranges  22®  to  24.5®  where  the 
back  branch  of  the  first  triplication  causes  some  clear  features  as  indi¬ 
cated  by  the  arrows.  Records  from  both  the  European  events  and  the  Idaho 
event  are  fit  well  by  the  model.  Figure  14  compares  data  to  synthetics 
for  A  ■  24.5®  to  28®.  The  model  fits  all  of  the  observations  from  the  U.S. 
and  Europe  reasonably  well.  The  records  of  the  4  new  events  being  studied 
did  not  contain  clear  features  associated  with  the  back  branch  of  the 
second  triplication.  For  this  reason  two  records  from  the  Borrego  Mountain 
earthquake  which  do  contain  such  an  arrival  were  modeled.  The  records  are 
marked  "BO"  in  the  figure  and  the  arrival  is  indicated  by  the  arrows.  This 
arrival  was  originally  modeled  with  T7  by  Burdick  and  Helmberger  11978], 
but  T9  fits  it  just  as  well.  Thus  we  have  included  in  the  data  set  clear 
features  associated  with  the  forward  and  back  branches  of  both  tripli¬ 
cations  and  have  modeled  them  with  models  S8  and  T9.  No  strong  evidence 
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has  been  found  that  the  sizes  of  the  discontinuities  or  the  basic  shape  of 
the  triplication  curves  vary  systematically  in  North  America  and  Europe. 


RESULTS  AND  DISCUSSION 


The  process  of  searching  for  large  scale  lateral  variations  in  upper 
mantle  structure  using  the  synthetic  seismogram  method  has  only  begun. 
Many  continental  as  well  as  oceanic  areas  have  not  yet.  been  studied  by  any 
method.  However,  there  are  some  important  results  which  are  true  for  this 
particular  comparison.  There  is  no  consistent  evidence  for  broad  scale 
differences  in  such  features  as  the  size  of  the  major  discontinuities, 
their  relative  locations  or  the  strength  of  the  velocity  gradients  between 
them.  There  is  no  evidence  for  additional  discontinuities.  In  the 
quintuplication  region,  where  the  effect  of  the  upper  mantle  discontin¬ 
uities  is  strongest,  we  find  very  similar  upper  mantle  delta  function 
responses  for  North  America  and  Europe  (Figure  7).  Most  of  the  evidence 
for  lateral  variations  in  upper  mantle  structure  can  be  explained  in  terms 
of  differences  above  250  km. 

Given  and  Helmberger  119801  completed  a  waveform  analysis  study  for 
northwestern  Eurasia  which  gave  results  very  similar  to  the  ones  presented 
here.  Figure  15  compares  their  preferred  model  K.8  with  S8  of  this  study. 
The  sizes  of  the  discontinuities,  their  relative  separation  and  the 
strengths  of  the  gradients  between  them  are  virtually  the  same.  The  inset 
shows  that  the  shape  of  the  triplications  in  Che  travel  time  curve  are 
nearly  the  same  for  the  two  models.  However,  there  is  a  difference  in  the 
absolute  depth  to  the  discontinuities  of  about  20  km  causing  a  relative 
shift  in  the  triplications  of  about  l®.  The  differences  in  depth  to  the 
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discontinuities  may  not  be  significant.  Different  structures  were  adopted 
in  the  two  studies  for  the  crust  and  lid.  The  apparent  variations  in 
absolute  depth  could  be  partially  or  completely  absorbed  in  the  top 
portions  of  the  models.  The  shift  in  the  distance  range  of  the  tripli¬ 
cation  by  1°,  on  the  other  hand,  is  grounded  directly  in  the  observations. 
It  also  corresponds  with  the  types  of  lateral  variations  found  -in  this 
study.  There  is  a  definite  possibility  that  systematic  variations  in 
absolute  depth  to  the  transition  zones  do  occur  at  about  the  level  indi¬ 
cated  in  Figure  13.  The  studies  of  Jordan  [1979]  may  have  resolved  this 
type  of  variation.  If  the  absolute  depth  to  the  major  discontinuities  does 
change  on  a  broad  scale,  but  their  relative  shape  and  position  remains 
constant,  then  vertical  delay  times  would  change,  but  the  shapes  of  the 
triplication  curves  would  not.  This  type  of  lateral  variation  model  could 
reconcile  all  the  observations,  although  it  would  imply  lateral  density 
gradients.  The  differences  in  the  depths  to  the  discontinuities  would  also 
have  to  be  larger  than  those  in  Figure  15  to  explain  the  magnitude  of  the  S 
wave  residual  differences  between  continents  and  oceans.  A  final  impor¬ 
tant  point  is  that  the  S  wave  vertical  travel  time  residuals  cited  by 
Jordan  [19791  varied  by  only  about  1.5  seconds  between  stable  and  active 
continent.  In  P  waves,  this  would  correspond  to  a  difference  of  only  about 
1  second.  Models  S8  and  T9  can  account  for  this  difference. 

One  possible  shortcoming  of  the  upper  mantle  models  for  northwestern 
Eurasia  of  both  Given  and  Helmberger  [1980]  and  King  and  Calcagnile  [1976] 
is  that  they  are  based  on  travel  time  data  that  may  be  inaccurate.  Both 
models  were  constrained  to  fit  ISC  travel  times  from  Soviet  nuclear  tests. 
These  observations  dictate  that  travel  times  from  25°  to  30*  should  be 
very  close  to  JB.  Rodean  [1979]  points,  however,  that  the  estimated  origin 
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times  for  these  events  are  always  2  or  3  seconds  before  the  minute.  He 
concludes  that  the  distribution  of  stations  around  the  USSR  combined  with 
the  use  of  JB  travel  time  curves  by  ISC  in  locating  events  has  caused  this 
systematic  error  in  origin  time.  Many  studies  show  that  stable  continental 
travel  times  are  faster  than  JB  from  25®  to  30°  [Hales,  1972;  Carder  e£ 
al .  ,  19761,  and  S8  has  been  designed  to  fit  this  criterion.  Figure  13 
shows  that  if  model  K8  was  adjusted  so  that  it  was  2  seconds  faster  for 
rays  that  bottom  below  800  km,  it  would  begin  to  resemble  S8  even  more. 

The  previous  work  of  King  and  Calcagnile  [  1976]  and  England  et  al . 
11977,  19781  suggested  that  major  lateral  variations  in  the  sizes  of  the 
upper  mantle  discontinuities  existed  in  the  continental  regions  sur¬ 
rounding  NORSAR.  Waveform  analysis  studies  hax'e  failed  to  substantiate 
this  result.  The  King  and  Calcagnile  model,  KCA,  contained  an  unusual  low 
gradient  zone  above  the  first  discontinuity  and  a  relatively  large  (7£) 
jump  in  velocity  for  that  transition.  The  work  of  Given  and  Helmberger 
[19801  has  shown  that  these  unusual  features  are  not  all  representative  of 
the  mantle  structure  of  northwestern  Eurasia.  The  synthetics  from  KCA  do 
not  even  satisfy  waveform  observations  from  WWSSN  stations  very  close  to 
NORSAR.  The  low  gradient  zone  above  420  km  is  either  a  very  localized 
feature  or  an  artifact  caused  by  some  pecularity  of  the  shallow  structure 
beneath  NORSAR.  England  et_  al .  [  1977]  concluded  that  major  differences 
existed  between  Europe  and  northwestern  Eurasia  because  they  also  failed 
to  find  evidence  of  the  unusual  structure  in  KCA  above  420  km.  Since  this 
structure  represents  at  best  a  very  local  feature  and  since  very  similar 
models  have  been  found  for  the  two  regions  from  waveform  studies 
(Figure  15),  few  reasons  remain  to  suspect  deep  lateral  variations  in 


structure  exist  between  them.  It  would  appear  from  these  comparisons  at 
least  that  the  relative  shape  of  the  velocity  profile  through  the  trans¬ 
ition  region  is  approximately  as  shown  in  Figure  15  beneath  North  America, 
Europe  and  northwestern  Eurasia.  The  predominant  type  of  lateral  vari¬ 
ation  in  wave  propagation  is  a  shifting  of  the  triplications  in  distance 
caused  by  laterally  varying  structure  above  250  km  (Figure  8)  or  variation 
in  the  depth  to  the  transitions  (Figure  15). 
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TABLE  2 


Velocity  Models 


Depth  to  Bottom  Depth  to  Bottom 

of  Layer  (km)  a  (km/sec)  of  Layer  (km)  a  (km/sec) 


T9 


S8 


20.0 

6,00 

30.0 

6.  50 

35.0 

6.50 

50.0 

7.10 

45.0 

7.95 

60,0 

8.  10 

55.0 

8.00 

70,0 

8.20 

65.0 

8.05 

80,0 

8.24 

73.  1 

7. 75 

90,0 

8,24 

89,  8 

7.72 

100,0 

8,25 

103.3 

7.70 

120.  0 

3.25 

117.5 

7,73 

130.  0 

8-26 

127.0 

7,77 

140,0 

8.  26 

135.7 

7.83 

150.0 

8.  27 

144,4 

7.88 

160.  0 

8.  28 

150.5 

7,94 

170  0 

8.  29 

156,4 

8-03 

180.0 

8,  10 

162.3 

8.11 

i90,0 

8.  32 

168.2 

8.19 

200,  0 

8,  34 

180.0 

8-25 

210.0 

8.  36 

190,0 

8,28 

220  0 

8.  18 

200,0 

8.30 

230.0 

8.  40 

210.0 

8,33 

240.0 

8  42 

220.0 

8.35 

250.0 

8.44 

230.  0 

8.  38 

240.0 

3.41 

250.0 

8.43 

Both  T9  and  S8 

260.0 

8,46 

270.0 

8.48 

280.0 

8.51 

290.0 

8.  54 

300.0 

8,  56 

310.0 

8.  59 

320.0 

8.61 

330.0 

8.64 

340.0 

8.67 

350.0 

8.69 
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Depth  to  Bottom  Depth  to  Bottom 

of  Layer  (km)  a  (km/sec)  of  Layer  (km)  a  (km/sec) 


Both  T9  and  S8 


360.0 

8.72 

370.0 

8.74 

380.0 

8.77 

386.0 

8.79 

395.0 

8.81 

400.0 

9.24 

410.0 

9.27 

420.0 

9.31 

430.0 

9.35 

440.0 

9.39 

450.0 

9.43 

460.0 

9.46 

470.0 

9.50 

430.0 

9.54 

490.0 

9.58 

500.0 

9.62 

510.0 

9.65 

520.0 

9.69 

530.0 

9.73 

540.0 

9.77 

550.0 

9.80 

560.0 

9.84 

570.0 

9.88 

580.0 

9.92 

590.0 

9.96 

600.0 

10.00 

610.0 

10.03 

620.0 

10.07 

630.0 

10.11 

640.0 

10.18 

650.0 

10.25 

655.0 

10.  76 

665.0 

10.78 

675.0 

10.80 

685.0 

10.82 

695.0 

10.84 

705.0 

10.86 

715.0 

10.88 

725.0 

10.90 

739.0 

10.92 

745.0 

10.94 

755.0 

10.96 

765.0 

10.98 

775.0 

11.00 

785.0 

11.02 

795.0 

11.04 

Both  T9  and  S8 


805.0 

11.06 

815.0 

11.08 

825.0 

12.00 

835.0 

11.12 

845.0 

11.14 

855.0 

11.16 

865.0 

11.18 

885.0 

11.22 

905.0 

11.26 
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FIGUPE  CAPTIONS 

The  locations  of  the  seismic  events  (squares)  used  in  the 
study  of  North  America  are  shown  with  respect  to  the  WWSSN 
and  Canadian  stations  (circles).  The  shaded  area  is  the 
previously  modeled  tectonic  part  of  the  continent. 

The  locations  of  the  two  seismic  events  (triangles)  used  in 
the  study  of  Europe  are  shown  with  respect  to  the  WWSSN 
stations  (circles). 

Observed  records  (heavy  line)  from  the  Idaho  event  from  dis¬ 
tance  ranges  between  30®  and  90®  are  compared  to  the  corres¬ 
ponding  synthetics  (light  line).  The  data  and  synthetics 
shown  here  and  throughout  are  normalized  to  peak  amplitude. 

Observed  records  (heavy  line'  from  the  Bermuda  event  from 
distance  ranges  between  30®  and  90°  are  compared  to  the 
corresponding  synthetics  (light  line). 

Observed  records  (heavy  line)  from  the  first  European  event 
(El )  from  distance  ranges  between  30°  and  90®  are  compared 
to  the  corresponding  synthetics  (light  line'. 

Observed  records  (heavy  line)  from  the  first  European  event 
|'E1)  from  distance  ranges  between  30®  and  90®  are  compared 
'  :■  the  corresponding  synthetics  (light  line). 


The  two  layered  velocity  models  representing  stable  and 
tectonic  continental  structure  are  compared.  The  only  diff¬ 
erences  in  the  models  occur  above  250  km. 

The  triplication  curves  for  the  stable  and  tectonic  contin¬ 
ent  models  are  plotted  along  with  observed  data. 

The  delta  function  responses  for  models  S8  and  T9  are  com¬ 
pared  in  the  distance  ranges  where  the  differences  in  the 
low  velocity  zone  continue  to  have  a  large  effect.  The 
responses  have  been  smoothed  with  a  gaussian  filter  to 
remove  the  noise  due  to  model  layering. 

The  observed  seismograms  from  ranges  between  14°  and  18® 
are  shown  in  heavy  line.  To  the  left  is  shown  the  synthetic 
for  the  source  except  for  the  mantle  response  followed  by 
the  mantle  model  delta  function  response  and  then  the  final 
synthetic  (all  in  light  line)  beneath  the  observed  wave. 
Synthetics  are  shown  for  both  S8  and  T9.  The  arrows  indicate 
clear  features  associated  with  Che  triplication  arrival.  B 
indicates  records  from  the  Bermuda  earthquake,  I  from  Idaho, 
El  from  the  3/4/67  European  and  E2  from  the  11/30/67 
European. 

The  observed  seismograms  (heavy  line)  from  ranges  between 
19®  and  22®  are  compared  to  synthetics  (light  line).  The 
arrangement  of  data  and  synthetics  is  the  same  as  in 


Figure  10.  The  arrows  indicate  features  associated  with  the 
arrival  of  the  reflection  from  the  lower  discontinuity. 

Observed  P  waves  (heavy  line)  are  compared  to  synthetics 
(light  line)  computed  for  a  slightly  shifted  range.  The 
source  of  the  MAL  record  has  been  changed  but  the  distance  is 
correct.  The  arrangement  of  data  and  synthetics  is  the  same 
as  in  Figure  10. 

The  observed  seismograms  (light  line)  from  ranges  between 
22®  and  25®  are  compared  to  synthetics  (light  line).  The 
arrangement  of  data  and  synthetics  is  the  same  as  in 
Figure  10<  The  arrows  indicate  features  associated  with  the 
reflection  from  the  upper  discontinuity. 

The  observed  seismograms  (heavy  line)  are  compared  to  the 
synthetics  (light  line)  for  ranges  between  24®  and  29®.  The 
arrangement  of  data  and  synthetics  is  the  same  as  in 
Figure  10,  "BO”  indicates  records  from  the  Borrego  Mountain 
earthquake . 

The  layered  velocity  models  for  stable  continent  and  north¬ 
east  Eurasia  are  compared.  The  travel  time  curves  in  the 
inset  are  plotted  on  shifted  time  and  distance  scales. 
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THE  FOCAL  MECHANISMS  OF  THE  GAZLI,  USSR  EARTHQUAKES1 

2 

M.  J.  Kristy  ,  L.  J.  Burdick  and  D.  W.  Simpson 
Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  New  York  10964 

ABSTRACT 

Two  Mg  -  7  earthquakes  struck  in  April  and  May  of  1976  near  the  town  of 
Gazli,  Uzbekistan  in  the  Soviet  Union  causing  extensive  damage.  The  focal  mech¬ 
anisms  of  the  two  events  have  been  determined  by  modeling  the  long  period  body 
waves  with  synthetic  seismograms.  Both  events  have  large  components  of  thrust 

faulting,  but  they  vary  measurably  in  azimuthal  orientation.  The  moment  of  the 

26  26 
first  event  was  2.1  x  10  dyne-cm  and  of  the  second  was  2.4  x  10  dyne-cm.  The 

stress  drops  were  around  SO  bars  which  is  typical  for  intraplate  events. 

The  surrounding  region  had  been  assigned  a  low  seismic  risk  based  on  recent 

patterns  of  seismicity.  Historic  records  from  the  ancient  city  of  Bukhara  show 

that  earthquakes  in  the  magnitude  range  6  to  7  had  actually  occurred  in  the  region 

several  times  prior  to  the  turn  of  the  century.  The  Gazli  events  appear  to  have 

j 

occurred  on  a  pre-existing  zone  of  weakness  that  is  evidenced  by  the  geologic 


record. 
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INTRODUCTION 

Two  large  earthquakes  occurred  in  April  and  May  1976  in  the  Kyzyl-Kum 

desert  of  western  Uzbekistan,  USSR,  approximately  400  km  southwest  of  Tashkent 

(Figure  1).  The  first  event  of  2:40:23.9  April  8,  1976  was  located  at  40.3°N,  63.72°E 

(ISC)  about  30  km  northeast  of  the  small  town  of  Gazli.  It  had  a  magnitude 

M  =7.1  (PAS).  At  2:58:41.1  May  17,  1976,  the  second  event  occurred  20  km  west 

of  the  first  at  40.35°N,  63.45°£  with  a  magnitude  M  =7.1  (PAS).  Deaths 

s 

apparently  occurred  and  extensive  damage  was  reported  at  Gazli  and  in  nearby  gas 
fields  (Aptekman  et  al.,  1978;  Pletnev  et  ah,  1977). 

The  Gazli  earthquakes  have  several  unusual  aspects  that  are  of  special 
interest.  The  occurrence  of  two  large  earthquakes  of  the  same  magnitude  at 
essentially  the  same  location  within  six  weeks  of  one  another  does  not  fall  into  any 
of  the  usual  fores  hock-mains  hock-aftershock  patterns.  A  high  level  of  activity 
occurred  between  and  after  the  two  shocks,  but  the  active  area  was  small 
compared  to  earthquakes  of  similar  magnitude  (Aptekman  et  al.,  1978). 

These  earthquakes  occurred  in  an  intraplate  region  that  had  previously 
exhibited  relatively  low  seismic  activity.  Figure  2a  shows  all  of  the  events  from 
the  Kondorskaya  and  Shebalin  (1977)  catalogue  with  magnitudes  greater  than  7 
from  the  years  1900  to  the  present.  There  was  certainly  little  evidence  from  these 
events  that  large  earthquakes  would  occur  so  far  westward  in  the  Eurasian  block. 
Figure  2b  shows  the  current  pattern  of  smaller  seismicity  including  all  shallow 
(<100  km)  events  from  the  ISC  catalogue  from  1964  to  1977  that  were  located  with 
more  than  10  stations.  This  seismicity  again  gave  little  indication  of  the  potential 
of  the  Gazli  area.  Also  shown  in  2a  are  the  seismic  zoning  contours  that  had  been 
adopted  prior  to  the  earthquake.  The  area  was  assigned  a  low  seismic  risk 
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Medvedev,  1976)  with  a  maximum  expected  intensity  of  V  as  measured  on  the  1964 
MSK  scale  (roughly  equivalent  to  the  Modified  Mercalli  scale).  The  actual 
intensities  experienced  during  the  1976  earthquake  were  near  X,  and  an  accelero- 
graph  installed  before  the  May  17  event  recorded  a  maximum  vertical  acceleration 
cf  1.3  g  (Pletnev  et  al.,  1S77). 

Even  though  the  two  earthquakes  are  separated  by  only  20  km  (Aptekman  et 
al.,  1978),  we  present  data  below  showing  that  the  focal  mechanisms  are  substan¬ 
tially  different.  In  agreement  with  most  other  fault  plane  solutions  for  this  region 
(e.g.,  Shirkova,  1967;  Das  and  Filson,  1975;  Nowroozi,  1972;  Molnar  et  al.,  1973;  Ni, 
1978),  the  main  component  for  both  events  is  thrusting;  however  the  azimuth  of  the 
nodal  planes  differs  by  more  than  60°.  Much  of  the  broad,  diffuse  zone  of 
intraplate  seismicity  north  of  the  India-Eurasia  plate  boundary  can  be  interpreted 
as  occurring  along  pre-existing  zones  of  weakness  separating  stable  continental 
blocks  (e.g.,  Molnar  and  Tapponnier,  1975;  Das  and  Filson,  1975).  In  some  areas, 
these  zones  can  be  identified  based  on  geological  data  and  background  seismicity. 
Such  correlations  are  difficult  in  the  Gazli  area  because  of  the  Quaternary  cover  of 
the  Kyzyl  Kum  sands  and  of  the  paucity  of  recent  seismicity.  Historical  data  from 
the  ancient  city  of  Bukhara  do  show  that  large  earthquakes  have  occurred  in  the 
region  over  the  last  800  years.  We  later  show  evidence  that  this  activity  may  be 
associated  with  the  westward  extension  of  a  zone  of  Hercynian  deformation  which, 
farther  to  the  east  in  the  southern  Tien  Shan,  follows  one  of  the  most  seismic  zones 
in  Central  Asia. 

The  Gazli  earthquakes  are  located  near  a  site  of  active  gas  production.  This 
connection  raises  the  possibility  that  they  may  have  been  induced.  Confirmation  of 
this  must  await  further  details  on  the  extent  and  location  of  the  pumping.  If  the 
connection  is  real,  however,  the  Gazli  earthqukes  would  be  the  largest  known 
induced  events. 
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THE  SOURCE  MODELS 


The  fault  plane  solutions  for  the  two  main  shocks  were  determined  by 
matching  the  WWS SN  long  period  records  of  the  P  and  SH  waves  with  synthetic 
seismograms.  A  full  description  of  the  techniques  involved  is  given  by  Helmberger 
and  Burdick  (1979).  Some  recent  applications  of  the  methods  to  similar  problems 
include  the  work  of  Cipar  (1980),  Bache  et  ah  (1980)  and  Ebel  et  al.  (1978). 
Fundamentally,  a  long  period  synthetic  seismogram  is  computed  by  convolving 
together  a  series  of  linear  time  dependent  operators.  The  first  and  most  important 
represents  the  response  for  a  point  double  couple  source  of  given  orientation  (fault 
plane  solution)  buried  in  a  layered  half  space.  The  computations  make  use  of 
Haskell  propagator  matrices.  The  time  function  of  the  source  is  assumed  to  be  a 
simple  triangular  pulse  independent  of  azimuth  and  ray  parameter.  Other 
operators  that  are  convolved  in  represent  the  effects  of  attenuation  and  the 
recording  instrument.  The  records  used  are  selected  from  the  ranges  30°  to  90°. 
This  avoids  complications  caused  by  earth  structure  because  the  bottoming  point 
for  the  rays  is  in  the  smooth  lower  mantle.  There  are  no  discontinuities  to  cause 
triplications  or  other  features  in  the  travel  time  curve.  The  records  modeled  were 
also  selected  on  the  basis  of  their  signal  to  noise  ratio  and  their  azimuthal 
coverage  of  the  source.  Perturbations  of  the  source  orientation,  the  time  function 
and  the  depth  in  the  half  space  are  made  by  trial  and  error  until  the  synthetic 
seismograms  provide  a  satisfactory  fit  to  the  observed  data.  The  observed  long 

I 

period  P  waves  are  simply  digitized  and  windowed  at  30  seconds.  The  observed  S 
waves  on  the  other  hand,  must  be  rotated  into  SH  and  SV  components.  This  is  a 
difficult  and  unstable  process.  SH  waves  which  were  small  with  respect  to  the 
corresponding  SV  signal  were  not  weighted  heavily  in  the  trial  and  error  fitting. 


Figures  3  and  4  show  the  P  and  SH  waveforms  respectively  for  the  April  8 
earthquake.  Figures  5  and  6  show  those  for  the  May  17  event.  The  April  8  event 
was  15  km  deep  and  was  modeled  with  a  triangular  source  pulse  with  rise  time 
2  seconds  and  fall  off  time  3  seconds.  The  May  17  event  was  also  15  km  deep  and 
the  source  pulse  model  had  a  rise  time  of  3  seconds  and  a  fall  time  of  4  seconds. 
Both  events  were  thrust  with  a  small  but  resolvable  component  of  strike  slip.  The 
azimuthal  orientation  of  the  fault  planes,  however,  was  definitely  different  for  the 
two  events.  The  P  waves  in  Figure  3  are  fit  well  by  the  model  but  they  do  not  vary 
strongly  with  azimuth.  The  P  wave  data  constrain  the  mechanism  to  be  a  dipping 
reverse  event  in  that  no  focal  planes  pass  through  the  center  of  the  sphere.  The 
strikes  of  the  planes,  however,  are  not  well  constrained.  Neither  the  data  or 
synthetics  vary  strongly  with  azimuth,  so  the  theoretical  solution  can  be  rotated 
about  a  vertical  axis  without  seriously  deteriorating  the  fit.  As  is  usual  for  dip  slip 
events,  the  SH  data  provide  a  much  stronger  constraint  on  the  strike  of  the  fault 
plane.  Figure  4  shows  the  S  data  displayed  around  the  lower  focal  hemisphere  for 
SH.  The  nodal  lines  for  this  case  are  not  planes  intersecting  the  sphere,  but  they 
do  divide  the  focal  sphere  into  positive  and  negative  sections.  The  equations  for 
the  SH  nodes  are  given  in  Burdick  (1977).  As  the  figure  shows,  there  must  be  a 
nodal  line  between  stations  PTO  and  BUL  and  also  between  GDH  and  MAT.  This 
constrains  the  azimuthal  orientation  of  the  solution.  The  precise  values  for  the 
fault  strike  and  the  rake  were  chosen  to  give  the  best  match  to  the',  details  of  the 

j 

waveform.  For  the  P  waves  this  amounts  to  matching  the  relative  sizes  of  the  P, 
pP  and  sP  phases  as  they  vary  with  azimuth.  For  the  S  waves  the  relative  sizes  of 
S  and  sS  are  matched. 

The  fault  planes  for  the  May  17  event  are  constrained  by  the  data  in  much 
the  same  way.  The  P  waves  show  the  mechanism  to  be  a  dipping  thrust  and  the  SH 
waves  constrain  the  strike.  Since  it  may  be  surprising  to  some  that  the  two  shocks 
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had  substantially  different  mechanisms,  the  evidence  that  shows  this  is  worth 
covering  carefully.  In  Figure  3,  the  P  wave  observations  to  the  east  (SEO,  SHK, 
BAG,  MAN,  DAV)  show  a  relatively  strong  first  motion  upward  caused  by  a  strong 
direct  P  phase.  The  corresponding  stations  in  Figure  5  show  a  nodal  first  motion, 
particularly  MAT,  SEO,  HKC,  DAV  and  LEM.  In  the  fault  plane  solution  in  Figure  6 
these  stations  are  very  close  to  the  P  node,  so  the  synthetics  match  the  effect 
quite  well.  The  evidence  from  the  SH  waves  that  the  two  focal  mechanisms  were 
different  is  even  stronger.  The  waveforms  at  KTG,  GDH  and  PRE  in  Figure  4  are 
much  different  than  the  corresponding  ones  at  KTG,  COL  and  PRE  in  Figure  6. 
The  comparison  between  GDH  and  COL  is  not  strictly  as  valid  as  the  other 
comparisons,  but  we  argue  that  since  they  are  so  close  in  location  on  the  focal 
sphere  their  waveforms  should  be  very  similar. 

Aptekman  et  al.  (1978)  and  Bezrodny  (1979)  present  fault  plane  solutions  for 
the  Gazli  events  based  on  first  motion  data  from  teleseismic  stations  and  local 
Soviet  networks.  Although  there  are  major  differences  between  the  two  papers 
they  agree  that  the  two  Gazli  events  had  substantially  different  mechanisms. 
Their  focal  mechanisms  are  consistent  for  the  most  part  with  the  results  of  this 
investigation  with  the  exception  of  Bezrodny's  solution  for  the  April  event. 
However,  his  result  is  very  inconsistent  with  the  SH  waveform  data.  A  very 
thorough  investigation  of  the  mechanism  of  the  May  17  event  has  recently  been 
completed  by  Hartzell  (1980).  His  fault  plane  solution  is  similar  to  the  one  in 
Figure  5  but  rotated  clockwise  by  15°.  He  found  that  the  slight  rotation  was 
required  by  the  surface  wave  data.  He  used  a  finite  fault  model  to  further  improve 
the  fit.  A  careful  comparison  of  his  finite  fault  synthetics  with  those  in  Figure  5 
indicates  how  difficult  it  is  to  resolve  the  effects  of  source  finiteness  from  long 
period  data.  Pearce  (1980)  has  modeled  one  of  the  aftershocks  of  the  April  event. 
He  found  that  the  mechanism  was  a  near  45°  dip  slip  reverse  fault  with  a  strike  slip 
component. 
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The  depths  of  the  events  were  constrained  by  matching  the  relative  timing  of 
the  depth  phases.  A  good  discussion  of  the  relative  sizes  and  polarities  of  the 
depth  phases  for  dipping  thrust  events  can  be  found  in  Langston  and  Helmberger 
(1975).  3riefly,  in  all  the  P  waves  modeled  here,  the  direct  P  is  a  small  positive 
arrival  followed  by  strong  pP  and  sP  arrivals  which  are  negative.  This  causes  the 
small  first  swing  -  large  second  swing  shape  of  the  P  waves  in  Figures  3  and  5.  The 
interaction  of  the  S  and  sS  phases  is  much  more  variable,  but  the  KTG  record  in 
Figure  6  shows  a  case  where  the  S  arrival  is  just  a  small  negative  arrival  followed 
by  a  strong  positive  sS  arrival.  The  relative  timing  of  the  phases  is  well  matched  in 
the  synthetics.  The  actual  values  of  depth  obtained  depend  of  course  on  the  crustal 
model  used.  Some  information  on  the  velocity  structure  near  Gazli  is  available 
from  the  work  of  Rodriguez  (1969)  and  Zunnonov  et  al.  (1974).  We  adopted  the 
model  of  Kosminskaya  et  al.  (1969)  since  it  seems  fairly  representative  (Table  I). 

The  moments  of  the  events  were  determined  from  the  observed  P  wave 

amplitudes  at  each  station  (Langston  and  Helmberger,  1975).  Figures  4  and  6  show 

the  moment  value  determined  from  each  station  next  to  the  corresponding 

26 

waveforms.  The  observed  values  for  the  April  8  event  range  from  0.7  x  10  to 

Off 

3.9  x  10  dyne-cm.  The  average  value  is  2.1  x  10  dyne-cm  with  a  standard 
deviation  of  0.9  dyne-cm.  The  average  observed  value  of  seismic  moment  for 
the  May  17  event  was  2.4  dyne-cm  with  a  range  of  0.4  to  4.8  x  it*'*"  dyne-cm  and  a 
standard  deviation  of  0.9  dyne-cm. 

Kanamori  and  Anderson  (1975)  show  that  interplate  and  intraplate  earth¬ 
quakes  segregate  on  a  moment-magnitude  diagram.  The  Gazli  events  fall  in  the 
center  of  the  distribution  of  intraplate  events.  Assuming  that  the  corresponding 
interpretation  of  stress  drop  from  Kanamori  and  Anderson  (1975)  is  correct,  we 
find  stress  drop  values  of  50  bars  for  both  events.  It  is  important  to  note,  however, 
that  the  relationships  between  stress  drop  and  other  fault  parameters  are  still  a 
point  of  controversy.  A  full  summary  of  the  source  models  is  given  in  Table  2. 
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The  largest  uncertainties  in  the  source  model  are  in  the  time  function  since  it 

trades  off  directly  with  the  Q  operator.  We  have  assumed  a  Futterman  (1962) 

operator  with  a  t*  of  1.0  s  and  a  t*  of  4.0  s.  We  could  have  used  different 
a  p 

values  and  absorbed  the  difference  directly  into  the  source  function.  Minster 
(1978)  has  pointed  out  that,  for  realistic  rheologies,  Q  must  have  some  frequency 
dependence.  We  find  that  our  synthetics  are  virtually  unchanged  if  we  use  his 
standard  linear  solid  Q  operator  with  tm  of  .2  s,  t*o  =  1.0  s,  and  t*  s  4.0  s. 
These  values  are  consistent  with  the  observations  of  Sipkin  and  Jordan  (1979).  If 
we  assume  the  above  values  for  t*  we  find  that  the  time  function  duration  is 
resolved  to  within  about  50%  and  the  values  of  rise  time  versus  fall  time  trade  off 
with  the  source  depth. 


TECTONICS  AND  SEISMICITY 


The  Gazli  earthquakes  are  well  removed  from  any  major  plate  boundary 
(Figure  1).  The  Arabian  and  Indian  plates  are  converging  on  the  Eurasian  plate,  but 
almost  1000  km  to  the  south  of  Gazli.  More  sophisticated  interpretations  of  the 
tectonics  involving  several  additional  plates  were  presented  by  Morgan  (1968), 
Le  Pichon  (1968)  and  Nowroozi  (1972)  among  others.  A  consistent  feature  in  all  of 
these  models  is  that  the  convergence  is  basically  north-south  and  that  the  Gazli 
region  is  not  associated  with  a  plateTooundary. 

Much  of  southern  and  eastern  Asia  shows  a  high  level  of  intraplate  seismicity, 
and  several  theories  have  been  proposed  to  explain  this  phenomenon.  Molnar  and 
Tapponnier  (1975)  compared  the  sense  and  orientation  of  faulting  in  Asia  to  the  slip 
lines  for  a  rigid  identor  in  a  plastic  medium  and  found  many  similarities.  Molnar  et 
aL  (1973)  and  Das  and  Filson  (1975)  suggest  that  zones  of  weakness  separate  Asia 
into  blocks  with  relatively  stable  interiors.  Sykes  (1978)  observed  that  in  many 
intraplate  areas  seismicity  tends  to  concentrate  along  zones  of  pre-existing 
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weakness  and  suggested  that  a  similar  phenomenon  is  occurring  in  Asia.  This 
concept  is  at  least  applicable  in  Soviet  Central  Asia  were  Nersesov  (personal 
communication)  has  found  the  seismicity  can  be  best  explained  in  terms  of  the 
movement  of  a  few  blocks. 

The  seismicity  and  major  faults  of  southwest  Asia  are  shown  in  Figure  2. 
Most  of  the  faults  shown  in  2a  are  evident  in  the  seismicity  pattern.  These  include 
the  Kopet  Dag  fault,  the  southern  Tien  Shan  fault  system,  the  Sharud  fault,  and  the 
Herat  fault  (Yanskin,  1966).  The  seismicity  between  the  Herat  and  southern  Tien 
Shan  faults  is  associated  with  the  Hindu  Kush.  The  segment  of  the  Talas-Fergana 
fault,  shown  in  Figure  2a  does  not  appear  to  be  seismically  active  at  present,  but 
an  extension  of  the  fault  to  the  southeast  is  active,  with  one  known  earthquake 
greater  than  magnitude  7  on  it. 

The  most  active  zones  of  seismicity  in  Figure  2b  follow  the  northern 
boundary  of  the  Alpine  orogenic  belt  (Figure  7)  through  the  Caspian  Sea,  along  the 
Kopet  Dag  and  Herat  faults  and  north  of  the  Pamirs.  In  many  parts  of  the  Tien 
Shan,  concentrated  zones  of  seismicity  can  be  related  to  edges  of  zones  of 
Hercynian  deformation,  e.g.,  along  the  Gissar-Kokshal  fault  in  the  southern  Tien 
Shan  (Gubin,  1967)  and  along  the  Talas-Fergana  and  main  Tien  Shan  faults  in  the 
northern  Tien  Shan  (Simpson  et  ai.,  1980).  Scattered  exposures  of  Hercynian  age 
continue  in  outcrops  northwest  of  the  southern  Tien  Shan  but  become  increasingly 
covered  by  Neogene  and  Quaternary  sands  (Markovskii,  1975;  Yanskin,  1966). 
Major  structural  features,  however,  extend  northwestward  through  Samarkand  and 
Bukhara.  Near  Gazli  these  split  into  two  trends,  one  continuing  along  the  same 
azimuth  and  passing  through  the  northern  Caspian  Sea  (Figure  7);  the  other  turning 
northward  to  join  the  Hercynian  structures  of  the  Ural  Mountains,  north  of  the  Aral 
Sea  (Beloussov,  1976). 
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The  Gazli  earthquakes  occurred  along  the  southern  boundary  of  the  Kyzyl- 
Kum  uplift  (Lukina,  1978).  Ulomov  (1974,  1979)  has  stressed  the  importance  of 
structural  trends  in  controlling  the  seismicity  of  the  Kyzul  Kum.  In  1974,  based  on 
surface  deformation  near  Tamdu-Bulak  (100  km  northeast  of  Gazli)  and  temporal 
changes  in  seismicity  in  western  Uzbekistan,  he  had  identified  the  Kyzul  Kum  area 
as  the  likely  site  for  a  large  earthquake  in  the  near  future  (Ulomov,  1974).  From 
'  >  geophysical  data,  several  faults  with  east-west  to  northeast-southwest  orientation 

have  been  proposed  for  this  region  (Lukina,  1978;  Hamrabaev  et  al.,  1976;  Zunnonov 
et  al.,  1974).  The  Kyzyl-Kum  uplift  provides  the  link  between  the  Hercynian 
%  deformation  in  the  Urals  and  the  Hercynian  deformation  in  the  southern  Tien  Shan 

4 

i  (Burtman,  1975;  Beloussov,  1976;  Kravchenko  et  al.,  1977).  As  observed  in 

i  "  1  " 

* 

]  Figure  7,  Hercynian  structures  in  the  southern  Tien  Shan  have  an  east-west 

*  orientation  and  the  Ural  Hercynian  belt,  north  of  the  Aral  Sea,  is  oriented  along  a 

north-south  trend.  The  connection  of  these  two  sections  of  the  Hercynian  orogeny 
has  subsided  and  is  an  arcuate  structure  that  continues  west  from  the  southern  Tien 
Shan,  northwest  along  the  Kyzyl-Kum  uplift  and  along  the  Aral  Sea  to  the  Urals 
(Beloussov,  1976).  We  propose,  therefore,  that  the  Gazli  earthquakes  occurred  on 
this  pre-existing  zone  of  weakness  that  follows  the  belt  of  Hercynian  deformation 
joining  the  Urals  and  Tien  Shan.  In  addition,  this  zone  is  especially  complex  near 
Gazli  where  it  undergoes  a  major  change  in  azimuth. 

i 

DISCUSSION 


The  1976  earthquakes  were  by  far  the  largest  events  in  the  area  in  recent 
times.  There  was  little  evidence  in  the  instrum entally  recorded  seismicity  nor 
obviously  in  the  more  extensive  data  considered  in  the  seismic  zoning  to  indicate  a 
high  seismic  potential  for  Gazli.  The  large  shocks  since  1900  had  all  been  located 


well  away  from  the  region.  The  few  earthquakes  located  near  Gazli  were  small  and 
occurred  prior  to  1960,  when  the  locations  were  relatively  inaccurate. 

A  somewhat  different  picture  emerges  if  we  consider  the  available  infor¬ 
mation  from  before  1900.  The  city  of  Bukhara  is  an  ancient  one  with  historical 
accounts  stretching  back  to  before  1200  A.D.  These  records  show  that  the  city  had 
been  struck  repeatedly  by  large  earthquakes.  Figure  8  shows  their  locations  as 
given  by  Kondoroskaya  and  Shebalin  (1977).  They  have  been  assigned  magnitudes 
from  6  to  7  based  primarily  on  intensity  data  from  Bukhara.  A  comparison  of 
Figures  2  and  8  shows  the  importance  of  the  historical  data.  There  is  a  clear 
indication  that  the  zone  of  seismicity  west  of  the  Tien  Shan  fault  system  reaches 
several  hundred  kilometers  farther  into  the  Eurasian  block  than  shown  by  the 
recent  data.  The  Bukhara  data  indicate  that  earthquakes  as  large  as  magnitude 
6  to  7  may  occur  along  this  zone  as  frequently  as  every  150  years.  Sykes  (1978) 
stresses  the  importance  of  identifying  such  zones  of  weakness  in  the  contiental 
crust  in  the  determination  of  seismic  risk. 

The  fact  that  the  two  Gazli  events  had  substantially  different  mechanisms  is 
well  resolved  by  the  data.  As  methods  of  determining  fault  mechanisms  have 
improved,  however,  this  type  of  behavior  in  multiple  shocks  has  been  found  quite 
commonly.  Some  particularly  analogous  cases  include  the  double  shock  thrust 
events  in  Friulli,  Italy  studied  by  Cipar  (1980),  the  sequence  north  of  the  Pamirs 

i 

studied  by  Jackson  et  aL  (1979)  and  the  1996  Karatag  earthquakes  in  Tadjikistan, 
west  of  Dushanbe  (Gubin,  1967).  Shebalin  (see  Aptekman  et  al.,  1978)  suggests  that 
such  double  shocks  are  characteristic  of  regions  in  which  tectonic  structures  do  not 
follow  simple  linear  trends.  In  fact,  following  the  April  8  Gazli  earthquake, 
Shebalin  had  identified  the  region  as  one  susceptible  to  multiple  shocks  and 
predicted,  based  on  comparisons  with  double  shocks  in  other  areas,  that  a  second 
large  earthquake  would  occur  during  the  second  half  of  May  (Aptekman  et  al.. 
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Sykes  (1978)  points  out  that  many  intraplate  earthquakes  occur  along  pre¬ 
existing  zones  of  weakness;  especially  where  these  zones  are  oriented  favorably 
with  respect  to  the  current  stress  direction.  In  regions  of  high  stress  where  no  such 
"favorably-oriented"  weak  zones  exist,  stress  may  be  relieved  in  pairs  of  earth¬ 
quakes  or  complex  single  events  rupturing  along  conjugate  zones,  oblique  to  the 
current  stress  direction.  We  note  that  the  nodal  planes  for  the  Gazli  earthquakes 
are  approximately  parallel  to  the  strike  of  the  two  trends  in  the  structural  fabric 
of  the  Gazli  area  (Figure  7).  Other  studies  (Shirokova,  1967;  Nowroozi,  1972; 
Molnar  et  al.,  1973;  Das  and  Filson,  1975)  show  a  general  northwest  compression 
along  the  USSR-Afghanistan  border  and  northeast  compression  between  Iran  and 
the  USSR  (Figure  2a).  The  difference  in  stress  directions  inferred  from  the  two 
Gazli  earthquakes  reflects  the  strong  control  on  the  fault  planes  imposed  by  local 
geological  structure.  Hence,  although  the  thrusting  mechanisms  are  indicative  of 
strong  horizontal  compression  and  are  in  agreement  with  previous  studies,  the 
inferred  directions  of  maximum  compression  from  Gazli  are  less  reliable  indicators 
of  regional  trends  than  the  more  consistent  directions  found  from  earthquakes  to 


the  south. 
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FIGURE  CAPTIONS 


Figure  1:  Major  plate  boundaries  and  their  relationship  to  the  Gazii  earth¬ 

quakes.  Closed  sawtooths  represent  ocean-continent  convergence 
and  the  open  teeth  are  continent-continent  convergence.  Single 
large  arrows  show  the  direction  of  relative  motion  between  two 
plates.  The  star  is  the  location  of  the  Gazii  earthquakes  and  the 
circle  is  the  triple  junction  for  the  Indian,  Arabian  and  Eurasian 
plates.  Plate  motions  and  boundaries  are  from  Jacob  and 
Quittmeyer  (1979). 

Figure  2:  Compressive  stress  orientations  and  the  relationship  between  the 

seismic  activity  and  seismic  zoning  in  the  USSR.  In  Figure  2a,  the 
solid  circles  are  instrum entally  located  earthquakes  with  magni¬ 
tude  >7  from  1900  to  1976.  The  heavy  solid  lines  are  major  faults: 
K  -  Kopet-Dag  fault,  ST  -  Southern  Tien  Shan  fault  system,  S  - 
Sharud  fault,  and  H  -  Herat  fault.  Thin  continuous  lines  show 
seismic  zoning  within  the  USSR  (Medvedev,  1976).  Bars  in  this 
figure  show  the  orientation  of  P-axes  from  published  focal  mechan¬ 
isms  (Shirokova,  1967}  Nowroozi,  1972)  and  from  the  Gazii  focal 
mechanisms  determined  in  this  study  (enclosed  in  the  box), 
'igure  2b  shows  the  same  area  and  includes  all  of  the  shallow 
(<100  km)  earthquakes  from  1964  through  March  1976  that  were 
located  by  the  ISC  with  10  or  more  stations. 
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Figure  3:  Fault  plane  solution  (lower  hemisphere)  and  corresponding  P  wave- 

forms  for  the  April  8,  1976  Gazli  earthquake.  Solid  circles  are 

compressions  and  open  circles  are  dilatations.  The  top  waveform  at 

each  station  is  the  observed  P  wave  and  the  lower  one  the  synthetic 

waveform.  In  the  lower  right  corner  of  each  box  is  the  body  wave 
26 

moment  {*10  dyne-cm)  for  each  station. 


Figure  4;  Fault  plane  solution  and  SH  waveforms  for  the  April  8,  1976 

< 

earthquake.  The  same  convention  is  used  as  in  Figure  3.  .  . 

Figure  5:  Fault  plane  solution  and  P  waveforms  for  the  May  17,  1976  earth¬ 

quake.  Notice  the  nodal  character  of  the  stations  MAT,  SEO,  HKC, 
and  DAV  in  comparison  to  similar  stations  for  the  April  8  earth¬ 
quake,  Figure  3. 


Figure  6;  Fault  plane  solution  and  SH  waveforms  for  the  May  17,  1976 
earthquake.  Notice  that  the  station  KTG  is  much  more  nodal  than 
that  observed  for  the  April  8  earthquake.  Figure  4. 


Figure  7:  The  relationship  between  the  Gazli  earthquakes  and  faults  as 

determined  from  geophysical  data  and  surface  outcrops  of 
Hercynian  age.  Thick  continuous  lines  are  "deep"  faults  and  lines 
with  dotted  interior  delineate  exposed  basement  of  Hercynian  age 
(Yanskin,  1966).  The  line  shading  represents  Alpine  deformation. 
Dashed  box  shows  the  area  of  the  inset.  Again  the  thick  lines  are 
faults  from  geophysical  data  (Zunnonov,  1974)  and  the  lines  with 
dashes  outline  local  outcrops  of  Hercynian  age. 
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Figure  8:  Historic  earthquakes  located  near  Bukhara  (hexagon)  from  intensity 

data  (Kondorskaya  and  Shebalin,  1978).  The  triangle  is  the  April  8, 
1976  event  and  the  square  the  May  17,  1976.  The  shaded  areas  are 
aftershock  zones  from  ISC  data.  The  solid  line  through  the  figure 
separates  the  region  according  to  the  seismic  zoning  (Medvedev, 
1976).  Roman  numerals  show  the  maximum  expected  intensity.  The 
historic  earthquakes  shown:  942  (M  =  6.7),  1390  (M  =  6.1),  1822 
(M=6.4),  1822  (M  =  7.0).  Vertical  and  horizontal  bars  represent 
inaccuracy  in  location  for  each:  earthquake.  ;  ^  ^  » 
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During  the  term  of  this  contract  v  >rk  on  Tectonics  of  Asia  (I’.I.  I.vnn  it.  Svl 
was  carried  out  in  the  following  areas: 

1)  A  study  of  the  focal  mechanisms  of  small  magnitude  earthquakes  in 
Central  Asia  using  surface  waves. 
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20.  ( con ' t . ) 


Work  on  Wave  f .  rm  Studios  of  the  Velocity  and  «}  Structure  of  the  'Unit  U-  (I’.I. 
lurry  d.  Burdick)  has  boon  carried  out  in  Lhe  following  areas: 

1)  Waveform  modeling  studies  of  the  structure  of  the  upper  mant le  have 
been  carried  out  in  both  the  U.S.  and  Europe.  The  results  are  compared  to 
those  of  other  regions  to  form  a  composite  picture  of  lateral  variation-,  in 
upper  mantle  structure. 

2)  A  study  of  the  seismicity  patterns  in  the  region  of  the  (!a;'li  earth¬ 
quakes  in  the  Soviet  Union  was  carried  out  to  attempt  to  understand  vhv  these 
unusual  earthquakes  occurred  where  they  did.  The  focal  mechanisms  of  the  two 
main  shocks  were  determined  hv  analysis  of  the  body  waveforms. 

1)  A  site  quality  survey  of  WWSSN  stations  in  Tndia  and  Southern  Asia 
has  been  made.  It  was  found  that  all  stations  near  the  continental  collision 
boundary  are  situated  on  dipping  structure  that  strongly  distorts  particle 
me-  *■  ion , 

4)  An  analysis  is  presented  of  the  importance  of  measurements  of  the 
relative  attenuation  of  1’  and  S  waves  to  determination  of  the  f requeue.'  de¬ 
pt  udance  of  0. 


_ L’ilCLAiaiLLiLLLD _ 

SECUUlTr  ci  ftr  ‘  *  *  ■  < 


